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Abstract 

In the dimer model, a configuration consists of a perfect matching of a fixed graph. If the underlying 
graph is planar and bipartite, such a configuration is associated to a height function. For appropriate 
"critical" (weighted) graphs, this height function is known to converge in the fine mesh limit to a Gaussian 
free field, following in particular Kenyon's work. 

In the present article, we study the asymptotics of smoothed and local field observables from the 
point of view of families of Cauchy-Riemann operators and their determinants. This allows in particular 
to obtain a functional invariance principle for the field; characterise completely the limiting field on 
toroidal graphs as a compactified free field; analyse electric correlators; and settle the Fisher-Stephenson 
conjecture on monomer correlators. 

The analysis is based on comparing the variation of determinants of families of (continuous) CR 
operators with that of their discrete (finite dimensional) approximants. This relies in turn on estimating 
precisely inverting kernels, in particular near singularities. In order to treat correlators of "singular" 
g local operators, elements of (multiplicatively) multi-valued discrete holomorphic functions are discussed. 
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1 Introduction 



The dimer model is a classical model of statistical mechanics; it consists in sampling uniformly (or with 
weights) perfect matchings of a bipartite graph. For planar graphs, Kasteleyn showed that the partition 
function can be expressed as the PfafRan of a properly signed (weighted) adjacency matrix for the graph. In 
the small mesh limit, this operator can be interpreted as a finite difference version of a 9-operator. In the 
"Tempcrleyan" case we will be considering, dimer configurations are in measure-preserving bijection with 
spanning trees on a related graph. 

Following Thurston, one may associate a height function to a dimer configuration on the square lattice 
(or another bipartite graph). As suggested by Benjamini, this height function can be understood in terms 
of windings of the associated spanning tree. 

In [26], Kenyon considers the scaling limit of the height function of dimers on the square lattice for 
planar domains with appropriate boundary conditions, and proves that in the scaling limit, it converges in 
distribution to a conformally invariant object for n-connected domains. This is based on precise asymptotics 
of the inverse Kasteleyn matrix. In [57], he identifies this limiting distribution as the classical massless 
Gaussian free field, in the simply connected case. 

In the present article, we are interested in several extensions of these results. The main point of view 
is that the invariance principles under consideration can be related in a natural way to Quillen's theory of 
families of Cauchy-Riemann operators ([35]). 

We outline the general strategy, first in the case of dimer tilings of the plane. Let h be the (discrete, 
random) height field and cj) a real compactly supported smooth test function. Then we shall see that 



E,exp,-./„Afl,.^ig|l,„., 



detiK) 

where if is a finite difference analogue of d, while Kg^ is an analogue of the Cauchy-Riemann operator 
d + (90), and (...) is a multiplicative correction which depends on 0. Here, h is piecewise constant on faces 
of the graph carrying the perfect matching. 

In order to take limits, it is convenient to consider Laplacian-type operators: 

det((ifgJ*ifaJ, 



EieM-^Jhm = (-) 
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In the continuum, there is a standard way to associate a "determinant" to such operators, namely 
regularisation, introduced by Ray and Singer f j39j). Let us be optimistic and assume for now that in 
the scaling limit, we can relate the above quantity to 

detc(^^) 
det(;(A) 

where A= d*d, Aa — {d + a)*{d + a). Then (a version of) the Quillen curvature formula ([31]) yields: 

du log(detcA,g^) = -^J m A {d<P) - / llV^/'f 
(since (j) is real). This is consistent with the defining property of the (normalised) massless free field: 

EpFieM^t / h^))^exp{ — / ||V0|n 

where cr is a scale constant. 

Now let us briefly discuss the situation for more complex topologies, in particular on a torus S. As 
shown by Kasteleyn, enumeration of dimer configurations by determinants is more complicate, involving 
spinor bundles. Also, the height function is only well-defined as a closed 1-form dh (alternatively, as an 
additively multiply- valued function h). The Hodge decomposition reads: 

dh — ujh + dh^ 

where oJh is a closed harmonic form (hence lying in a finite dimensional vector space ~ i?^(S); actually 
on a lattice A of this space due to model constraints) and ho is a single-valued function (defined up to an 
additive constant). We are interested in the limiting joint distribution of these objects. It is specified by 
the characteristic functional a i— >■ E(j3? J dhAa), where a ranges over (0, l)-forms (with associated Cauchy- 
Riemann operator d + a). The 9-Hodge decomposition gives a ^ Co + d4> where lj is an abelian form (a 
holomorphic 1-form). Now: 

dh Aa — J ujh Auj + J dho A d(j) 

Indeed, J ui A dho = - / {dui)hf = and / (Bcj)) Auh ^ ^ J {d(j>) A {uh + i * ujh) = 5 / {dc/)) A {uih + i * loh) = 
since dijjh — d * ujh — 0- The analysis will again be based on relating the family of CR operators a i-> 9 + a 
to its discrete counterpart, up to the level of {(,-) determinants. 

Beyond these "smooth" observables, for which convergence should quite robust, stricter conditions on the 
underlying graph (periodicity or isoradiality) allow to control "vertex" observables, such as electric correlators 
of type E(exp(i ajh{zj))) and monomer correlations (as in the Fisher-Stephenson conjecture). This will 
be addressed by identifying suitable families of CR operators and controlling the convergence of their finite 
difference approximations. 

New results obtained in this article include: functional convergence to the compactified free field for 
dimers on toroidal graphs; asymptotics of electric vertex correlators; and the Fisher-Stephenson conjecture 
for monomer correlations. In a follow-up article, extensions of these results, in particular to domains with 
boundary (bounded simply and multiply connected domains), will be considered. 

The article is organised as follows. Section 2 provides brief background on the continuous free field, the 
discrete dimer height field, and families of CR operators. In Section 3, we review the necessary results on 
isoradial dimers. Invariance principles in the plane are discussed in Section 4. The case of toroidal graphs 
is treated in Section 5. Section 6 describes a general surgery principle for (converging sequences) of discrete 
CR operators. Discrete multi-valued holomorphic functions and electric correlators are described in Section 
7. Monomers correlators are studied in Section 8. 
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2 Background and overview 



The main goal of this article is to express asymptotics of natural dimer observables in terms of a (com- 
pactified) free field; a key tool will be the study of families of (discrete and continuous) Cauchy-Riemann 
operators and their variational properties. In this section we provide a brief introduction to these objects. 

2.1 Free field 

The Gaussian (or massless) free field is a Gaussian field with covariance kernel given by the Laplacian Green's 
function (or a multiple thereof) . For a general introduction to the free field, see eg [101 [TOl ITl] . Here we will 
consider the free field on E, where E is the plane or a torus E = C/(Z + rZ). In both cases, the Laplacian 
has a non-trivial kernel consisting of constant functions. Consequently, the free field is only defined there up 
to an additive constant. 

Following Gross' abstract Wiener space approach, one may regard the free field as a random element 
of a Banach space, usually a Sobolev space modulo additive constants: _ff^^(E)/M (for some positive e). 
Note that elements of this space, and thus realisations of the free field, are distributions. In order to avoid 
quotienting by constant functions, one may consider its distributional total derivative, the current J = d(j). 

The (centered) free field is characterised by the fact that for any zero-mean test function ij] G C^(E), 
J^{tlj4>)dA is a centered Gaussian variable with variance g ip^^tpdA, where 17 > is a scaling constant 
(here A is the positive Laplacian). In other words, 

(exp(i / (V'(/))dv4) =exp(-| / ^pA-^^dA) 

where () denotes the expectation with respect to the free field. Instead of integrating against a smooth test 
function, it is often useful and natural to consider local "operators" , involving the local behaviour of the 
field at a certain number of selected points. For instance, one may consider current correlations, given by 
the Wick formula: 

n 

(J(zi) . . . J(z2n)) = ff" ^ JJrfz„^4^.Gs(2:Q.,Z^J 

{(Qi,ft),...,(a„,/J„)}i=l 

where the sum bears on the pairings {(ai, /3i), . . . , (q;„, /?«)} of indices in {1, . . . , 2n} (for an odd number 
of currents, the correlator vanishes due to the symmetry (f> -f-?- —4>). This expression makes sense when 
integrated by regular enough test functions ■(/'i(zi), . . . ,'02n(^2n)- 

After distributional derivatives of the field, the next local operators to consider are electric vertex oper- 
ators (|18[ll3j). which we write formally as exp{ia(j){z)) where the constant a is the charge. The difficulty is 
that, as (j) is not defined pointwise but as a distribution, it may not a priori be composed with the nonlinear 
function exp. A standard normalisation scheme for the correlator "(exp(^^- aj(l){zj))y^ consists in replacing 
exp(ia(/)(z)) with the (well-defined) exp{ia{TT5'^)^^ Jj^^^ (j)(z)dA{z)), letting 6 go to zero, and discarding 
the leading multiplicative term exp(— | J^j ctj (ttiS^)^^ J^^q Gci^i y)dA{x)dA{y)), which reflects the local 

blow-up of variance and not the geometry of the domain or the long range correlations. Alternatively, if 
is a field with continuous realisations and covariance kernel C, we have: 

(exp(i aj4>{zj)) = exp(-^ ^ ajakC{zj, Zk)) 

Here, the covariance kernel gGs blows up on the diagonal; to obtain a finite regularised correlator, we 
replace the coincident Green kernel "G'(a;, a;)" with the Robin kernel: G^{x, x) = r^{x) = lim^^a; G^{x, y) — 
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^ log \x — y\ (otherwise G^{x, y) = G^{x, y)) to obtain the expression: 

(: exp(i^aj0(zj) :) = exp(-| ^ Gs(2:i, Zfc)) 

where the colons recall the regularisation scheme, and we need J^j '^j = 0- In. the plane, this evaluates to 
the elementary expression: 



Zk 



j j<k 

Another set of "local operators" is given by magnetic correlators. These are disorder variables, ie they 
represent a modification of the state space rather than a modification of the state weights. Specifically, 
given marked points 21, . . . , 2„ and weights ai, . . . , q;„, one considers additively multivalued functions which 
increase by when cycling clockwise around ; one requires aj = so that the function is single- valued 
near infinity. This defines an affine state space with a single element of minimal (at least after normalisation) 
Dirichlet energy: 

00 = ^ aplog(z - zj) 

j 

In the Gaussian formalism, offsetting fields by a harmonic function (jjQ results in multiplying the partition 
function by exp(— /s |V</'op(^^) (see eg [HI for a discussion of free field partition functions). Here 
the Dirichlet energy is infinite, due to logarithmic singularities at the Zj's. A regularised Dirichlet energy 
is obtained by discarding the leading part of d(2 5)) l^'/'oP'^^- Representing by Oa{z) a magnetic 

charge a at z, we will consider the following regularised correlator: 



nO«,(^,):)=exp (^-^ 



1 p^Gg 



where 

rreg 



Again when E = C this correlator has a simple evaluation. Indeed, note that ipQ has the same Dirichlet 
energy as its (single valued) harmonic conjugate (/)o — ckj5Rlog(z — Zj). By Green's formula, 

/ |V^o|'dA = V / 4>odn4>od£ = -^y2{a,+0{S)){aj\og\5\+y2ak\og\zk-Zj\+0{6)) 

= -^C^(^])^og\6\ - ^^ajafclog|zj - Zk\ 



so that 



YlO„^iz,):)c^ll\z,-z,\ 

j 3<k 



Up to now, we have discussed the scalar free field, taking real values. An important f |181 113p variant 
is the compactified free field, taking values in the circle M/rZ; let us begin with an informal discussion. 
Given two Riemannian manifolds S and T, the classical harmonic mapping problem consists in finding a 
mapping (f> : S ^ T which minimises Dirichlet energy, for instance within a homotopy class. For example, 
the harmonic mappings C/L M/rZ are written z ^{zw), where ^{£iD) € rZ for all £ ^ L {L a lattice). 
In the quantised version of the problem, one consider mappings : S — >■ M/rZ (E a surface) with action 
functional given by the Dirichlet energy 

Si^) - ^ / \y^\'dVol 
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and distribution formally given by e~"''^'^^V<j). Observe that (j) can be decomposed as the sum of a harmonic 
mapping (j)^ and a scalar part (/)s : S — ^ M, projected on K/rZ; moreover the decomposition is unique up to 
additive constants. In terms of current, it amounts to considering the Hodge decomposition of J = d(j): 

J — dip — uJh + dho 

where is a closed harmonic form (the instanton component) and ho is a scalar field, well defined up to 
additive constant. Moreover, the action functional splits: 

^^(^J^u;hA*u;h + \Vho\^dVol^ 

It is thus natural to introduce a free field on E with compactification radius r as the data {oJh, ho) where the 
two components are independent. The 1-form uj^ has unnormalised distribution 

g-i f 

exp(-^ J UJh A *ujh)dfi{ujh) 

where /i is the counting measure on the lattice of harmonic 1-forms with periods integer multiple of r: 
UJh S rTj for any closed cycle A. The scalar component uj^ is a scalar free field with covariance kernel gG^ 
(and is defined modulo additive constant). 

Finally, let us mention that, similarly to Kramers- Wannier duality for the Ising model, the free field 
possesses an abelian duality which is a simple example of T- duality ([IHIIIS])- For a well-chosen normalisation 
of the field, T-duality inverts the compactification radius and exchanges electric (order) operators with 
magnetic (disorder) operators. 



2.2 Dimers 

For simplicity, we discuss here dimers on the square lattice Z^; however we shall subsequently consider the 
more general framework of isoradial graphs. For background on dimers, see |29| and references within. 

We will consider the square lattice itself, or some portion of it, or its quotient by a large scale lattice 
(yielding a graph embedded on a torus). A perfect matching or dimer configuration consists in a selection 
of edges such that each vertex abuts exactly one selected edge. For a finite graph, one may consider the 
uniform measure on such matchings. A Gibbs measure on matchings of the infinite volume graph 1? may 
be obtained as the weak limit of finite volume measures f |10l I31j ). The analysis of dimers, in particular in 
the fine mesh limit, relies extensively on Kasteleyn's Pfaffian enumeration of dimer configurations ([13]) for 
planar graphs. 

Following Thurston, to every dimer configuration on such a graph derived from 1? , one may associate an 
integer-valued (in some normalisation) height function h on its dual graph. Kenyon established in several 
set-ups that in the small mesh limit this height function converges to a free field ([IZ]); this has been 
extended in [TU |3T]. Specifically (for the square lattice), one can obtain a local central limit theorem, 
showing convergence of discrete current correlators to their free field limit. By integration against a test 
function, this gives the correct limit in law for the discrete height field, in some finite dimensional marginal 
sense. 

While the relationship between dimer height functions and the scalar free field was highlighted by these 
works, the relation with the compactified bee field becomes apparent when considering toroidal graphs. There, 
only the height current is well-defined, and the height function is additively multivalued (when tracing it 
along non-contractible cycles on the torus): as explained earlier, it splits into an instanton and a scalar 
component. In the case of the hexagonal lattice, the limiting distribution of the instanton marginal was 
worked out in [5 , in agreement with what is predicted by the compactified free field. This will be extended 
in the present article. 
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By analogy with vertex correlators for the free field, one may consider the asymptotics of discrete ob- 
servables of the type: (exp(z ajh(zj)), where (j) denotes the height function and the weights aj sum up to 
zero. The systematic study of these electric vertex correlators for dimer height functions has been initiated 
by Pinson (|37|) and will be further discussed here. 

A classical question for dimers, introduced by Fisher and Stephenson in [T7] , is the problem of monomer 
correlations. Specifically, the question is to estimate the variation of the partition function when a certain 
number of vertices are removed from the graph. A motivation for this question is a close analogy with Ising 
correlations, as illustrated by Hartwig's results ([H]) for monomer pair correlations on the diagonal. More 
recently, a series of articles by Ciucu (see [5] and references therein) has born on several variants of monomer 
correlations, in particular correlations of appropriate "islets" on the hexagonal lattice. In terms of height 
function, monomer insertions may be thought of as magnetic charges. Indeed, in the presence of monomers, 
the height function becomes additively multivalued, picking fixed additive constants when cycling around 
monomers (or larger defects). 

2.3 Families of Cauchy-Riemann operators 

We shall be mostly concerned with establishing convergence of dimer observables (height field, vertex cor- 
relators) to the corresponding (compactified) free field quantities. In all cases, the arguments will be based 
on the analysis of families of Cauchy-Riemann operators and their discrete counterparts. 

A Cauchy-Riemann operator in the sense of Quillen (|Mj) is an operator D : ^}P{L) — )■ ri'^'-'^(L) of type: 
D — dz{dg -\- a{z)) operating on sections of a complex line bundle L over a compact Riemann surface E, 
say. These operators constitute an affine (infinite dimensional) complex space. Given a metric on E,L, one 
can define an adjoint operator D* and the Laplacian-type operator D*D. The latter has a C-regularised 
determinant deti^{D* D) (eg [1]). In [3H], it is shown in particular how to evaluate the logarithmic variation 
of this determinant when D varies among Cauchy-Riemann operators. This is controlled by the asymptotic 
expansion of the inverting kernel of D near the diagonal. More precisely, if D is invertible, we have the near 
diagonal expansion (in the standard local coordinate): 

D-^{z, w) = D-^{z, w) + raiw) + 0{\z - w\) 

where is the parametrix: D~^{u,v) = exp{2i'^{a{w){u — v))) and a{w) = a{w)d'w {D~^ inverts 

the translation invariant operator that agrees with D at w). The curvature formula in |38| is based on the 
variational identity (specialised here to the simple case of line bundles on a flat torus): 

d f 

— logdet(;(L»t*A) J ratiw)dwAat 

where Dt = dz{d + at{z)) is a smooth parametric family of CR operators. This is itself a specialisation of 
the general variational formula for C determinants of Laplacian-type operators (0]). 

For instance, consider a torus S — C/(Z 4^ rZ). If we write a section s of a unitary line bundle as 
s(z) = /(z) exp(2i5ft(zA)), the d operator on this line bundle can be identified with the Cauchy-Riemann 
operator D = d + \dz, and A is a coordinate on the Jacobian of S. In this case, det(^{D* D) is the (square of) 
the Ray-Singer analytic torsion (P5]). Applying the Quillen variational formula in this case leads to Fay's 
evaluation of the Second Kronecker Limit formula (up to multiplicative constant). One may also consider 
CR operators written as D — d + ggdz; this corresponds to changing the metric on the trivial line bundle 
(ie scaling the fiber at z by e^^^^). In general, we are looking at d operating on a complex Hermitian line 
bundle and the variation of the associated C-determinant under a change of complex or Hermitian structure. 

In order to deal with vertex correlators, we will need to control more singular situations. For instance 
in the plane with n punctures S = C \ {zi, . . . , z„}, consider a fixed unitary character x of 7ri(E) and 
the associated complex line bundle L^^ , over E and d operator. There the punctures induce a logarithmic 
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divergence in the evaluation of det^ and an additional principal part regularisation is needed (see [55]). 
We will consider two types of variations: a moduli space (isomonodromic) variation (moving the punctures 
and keeping the character fixed) and a Jacobian variation (changing the character/complex structure of the 
line bundle with fixed punctures). 

On the discrete side, Kasteleyn ([53]) has shown how to evaluate the partition function of dimer config- 
urations on a planar graph as the determinant (in the bipartite case) of a nearest neighbour linear operator. 
This operator may be interpreted as a finite difference version of d, which is reflected in the asymptotic 
expansion of its inverse ( [261 1281 131j ) and underpins the analysis of the scaling limit of dimers. 

We will be relying systematically on this interpretation: in order to evaluate variations of the partition 
function, we will consider the associated modified Kasteleyn operator as a discrete version of a CR operator 
varying in a family. Variational formulae, both at the discrete and continuous level, involve the behaviour 
of the inverting kernel near the diagonal (after substracting the leading singularity, which is that of a 
translation-invariant operator). Convergence of observables will follow from convergence of these leading 
non-trivial coefficients in the short-range expansion of inverting kernels. For vertex correlators, this will 
involve a rather detailed study of discrete holomorphic functions with prescribed monodromy. 

The correspondence between observables and families of CR operators (or Hermitian line bundles) goes 
as follows: 

• instanton component o variation of the complex structure of the line bundle 

• scalar field variation of the Hermitian norm of the line bundle 

• variation of electric vertex correlators (: exp(i Q!j(^(zj) :) wrt position of insertions ^ family 
(zi,...,z„)^if^^) 

• variation of electric vertex correlators (: exp(j^^ aj(j){zj) :) wrt charges O family x ^(z ) 

Finally, for monomer correlations, for instance in the simplest case of pair correlations between points x 
and y, the relevant line bundle over C \ {x,y} has local sections s with monodromy (—1) around x,y and 

s{z) = O(y^fE^) near x,y. The variational analysis then consists in displacing one the punctures. 

In recent work ( 30 ), Kenyon employs related variational argument for rank 2 bundles in his study of 
the double dimer model. 

3 Dimers on isoradial graphs 

In this section, we are following the formalism of [35l [28l [7] for isoradial (critical) graphs. Notations and 
conventions are mostly as in [7j. 

3.1 Kasteleyn operator 

Consider a tiling A of the plane by rhombi with edge length 6. As faces have even degree, it is bipartite 
(ie there is a 2-colouring of vertices, say red and blue). Thus one may obtain two graphs from this tiling: 
the vertices of F are the blue vertices and are joined by an edge if they are vertices of the same rhombus; 
similarly the vertices F^ are the red vertices of A. The graphs F and F^ are dual. They are also isoradial 
in the sense that each face of F is inscribed in a circle of radius (5, the center of which is the corresponding 
vertex of the dual F^. The dual of A is denoted <0>. By abuse of terminology, we sometimes identify a graph 
with the set of its vertices. 
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We will work under the following assumption (see [7]): 



(4) : the lozenge angles belong to [?7o, tt — 770] for some fixed tjq £ (0, tt) 

Here 770 > is fixed once and for all; throughout, "absolute" constants may depend on 770 (but not on A etc). 

One can form a new bipartite graph M as follows: black (resp. white) vertices of M are the vertices 
(resp. centers) of rhombi in A. Edges of M are half-diagonals of the rhombi in A. We will be interested in 
perfect matchings of AI (or subgraphs of M). We refer to vertices of M as nodes. Black nodes corresponding 
to vertices of F (resp. vertices of F^) are vertex nodes (resp. face nodes); white nodes are edge nodes. 
Denote My ~ F, Mp ~ F^^, M]^ ~ <) the sets of vertex nodes, face nodes and edge nodes (ie white nodes) 
respectively, and Mb = My U Mp ~ A (black nodes). See Figure [l] 




Figure 1: (1) A portion of a rhombus tiling A. (2) Dual graphs F (solid) and F''' (dashed) obtained from A. 
(3) Corresponding bipartite graph M. 

Let Mb (resp. Mw) be the set of black (resp. white) vertices of M. The Kasteleyn operator K : C*^^ ^• 
£M„ gjygjj by (121]): 

(if/)(7i;) = ^if(7«,6)/(6) 

where K{w, bo) — §(^3 — 61) if bo, bi, 62, 63 are the black neighbours of w listed in counterclockwise order. 

For our purposes, it will be convenient to consider a real operator K : E*^^ — >• E*^"' which differs from 
K by the argument of its entries: K{w,b) — ±\K{'w,b)\, where the sign depends of a chosen orientation 
of {w,b) e Em, which now we describe (see [II|). Pick an arbitrary orientation of edges of F; then orient 
edges of F''' in such a way that if {xx') is an oriented edge of F, (yy') the dual oriented edge of , then 
{{xx'), (yy')) is a direct frame. Each edge of M is a half-edge of an edge of either F or F^ and inherits its 
orientation. Then K(it;,fe) = +\K{w,b)\ (resp. K(w,6) = -~\K{w,b)\) if [wb) is positively (resp. negatively) 
oriented. One readily checks that this yiekds a Kasteleyn orientation of A/, i.e. an orientation such that 
each face has an odd number of clockwise oriented edges on its boundary. 
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If w e Mw is on the oriented edge {xx') of T, set e'"'^'"^ = if^Efi ; also set e^"'^^'^ = 1 if & G T is a vertex 
node and e^''^''^ = i if 6 e FMs a face node. Then it is easy to check that 



for aU h, w, viz. K is obtained by composing K with diagonal operators. 

Translating results of [28] obtained for K in terms of K (see also fl}), we have: 

Theorem 1 ([2H])- There is a unique kernel K"^ : M*'^^ x M^^"' -> M smc/i </iot /or a// w G Miy, 

KK~^(., w) = Sio and K^^(6, w) = o(l) as |& — w| — > oo. 

2. If x,y,x',y' are black vertices in ccwise order around w, 

1 ^ ( y' 

K(?ii, a;)K~ (a;, w) = — arg 



27r \ y — X 



5. As \h -w\^ oo, 

/ iu(w) iv(b) 



7r(fo — w) J \\h — w 

Proof. See respectively Theorem 4.1, Theorem 4.2 and Corollary 7.4 in and also Theorem 2.14 in [7] 
for the general case. □ 

On M}^" , take the inner product given by counting measure and on R^^^ : 

{f,9)R'"w =^^J.o{w)f{w)g{w) 

w 

where /i^ (w) is the area of the rhombus of A containing w. With respect to these inner products, the adjoint 
operator K* ; M*^"' M*^^ has matrix elements: 

K*{b,w)=.pi<^{w)-^K{w,b) 

Then it is easy to check that (K*K)(6, b') = for aU 6 e T, 6' e T\ and if x ^ a;' in r and {yy') is the edge 
of r^' dual to {xx'), then 

^^^^^^"^'^'^ ^ (K*K)(x,a;) = ~J2 (K*K)(a;,x') 

(K*K)(y,y') = (K*K)(y,y) = - ^ (K*K)(y, y') 

so that with respect to the decomposition M^^^ M*^^ , K* K splits as Ar /Vt , where A is a graph 

Laplacian for appropriate edge weights (j28j). 

Define restriction operators Rb : C°iR^,C) and Rw ■ C^^'^^R^X) ^ K^"' by: 

(Rb0)(6) = 3?(e'''('VW) be Mb 

{Rwa){w) = 3fi(e-"^("')a/dz) «; e Mv^ 
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that restrict complex-valued functions (resp. (0, l)-forms) to black (rcsp. white) nodes. Note that these 
operators are only M-linear. If is a function on C, a direct computation shows that: 

(K(Rb<^))H = In (\y - y'Mw + u.^^) - ^{w - n.^^) + i\x - x'^w + m.^^) - ^{w - m.^^)' 



= 2fi^{w)md,(b{w)u) + o(l)) = 2^iM{f^w{^c|>){w) + o(l)) 

where {xx') is the oriented edge of F corresponding to w, {yy') the dual edge, u = e*''*^™^ and iJi<){w) is the 
area of the rhombus of A containing w (and thus of order 0(5^)). 

Let Xi,.. .,Xn be the neighbours of a vertex a; in F, in ccwise order (with cyclical indexing x„ — xa)- Let 
yk be the vertex of F''' corresponding to the face of F on the left-hand side of the oriented edge {xxk) and 
Wk S My/ the vertex corresponding to {xx^)- Then for a = Xdz, A G C constant, 

n n 

J2f^0{wk)e"'^'"''HRwa){wk) = ^iVk-i - 2/fc)5i(A(S^^)) 
fe=i fc=i 

= X! 4 (Kyk-i - yk){xk -x) + X{yk-i - yk){xk - x)) = ^'^IJ'oiwk) 

k=l k=l 

taking into account Xk — x = {yk-i ~ x) + {yk — x) and n<^{wk) = |(2/fc-i — yk){xk — x). Hence if tf) is 
continuous on a region Z) C C, 

^ ^,<,{w)e^^^^'^\Rw4>){^^) = \ I ^dA + o{l) = -- j aAdz + o{l) (3.1) 
since dA = dx A dy = ~^dz A dz. Similarly, 

n n ./ \ 

^Ar{x,Xk){xk - a;) = ^ ^ o'f ^^ _ -x) = Q 



fc=i 



^^{x,Xk){xk -xf = ^ ^^91^ ^- -xf = 

k=l k=l ^^^-'^ ^> 

which shows that the random walk on F associated to ZV converges to isotropic Brownian motion. 
We will also consider families of operators parameterized by a (0, l)-form a = a{z)dz: 

K„(u;, h) = K{w, h) exp ^23? j a j 

As before, if ^ is C^, we have 

( K„(i?B </)))( w) = ^^4^ URw{dcj)) + miauYRicj)) + 43fi(a(kt))3fi(i0) + o(l) 



2 

= 2^^(w) {Rwidcp + (pa) + o(l)) 

where = 4>{w), a = a{w), u = e^'^'^^\ 

In the case a = Xdz for some fixed A G C, we denote simply = K^df. We have Kx{w,b) = 
g-2K(A«;)K(i(;^5)e23*(A6) g^^^ consequently 

K^\b,w) t/' e-2^(^'')K-^(6,^i;)e2^(^™) 

satisfies K^K^^ = Id. 
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3.2 Height function 



Consider a finite bipartite graph S which is a subgraph of M bounded by a simple closed cycle on M . A 
perfect matching m of the bipartite graph S is a subset of edges of S such that every vertex of S is incident 
to exactly one edge in m. The weight of a matching m is the product of the weights of the edges present in 
m: 

w{m)= n \Hb,w)\ 

The partition function of the model is the sum of these weights over ah possible perfect matchings of M\ 

n 

iTi {bw)^m. 

We will consider the probability measure on perfect matchings m of S (if they exist) given by P(m = mo) = 
w{mo)/Z~^. 

In order to have at last one matching, it is necessary that ~ S H Afs and 'E.w = S H Mw have the 
same number of vertices. Let us also denote by K : E"^ M""' the restriction of K : M*^^ M*^"', ie 
(K/)H = EbeH«K(6,«;)/(&)- 

The fundamental result, due to Kasteleyn (|23|). is the following determinantal enumeration formula: 

Z = ±dct(K) 

(this is where the Kasteleyn orientation condition is required). Notice that K : M"^ — >■ M""' is not an 
endomorphism; however, Mr^ and M""' have canonical bases (up to permutation) with respect to which 
this determinant is evaluated (up to sign). Also, det{K) = det(K)n„g- e*''^"'' Ylba^ ^'''(^^ so that Z = 
\det{K)\. 

Moreover, we can also evaluate the partition function for arbitrary complex edge weights. If we replace 
the positive weight |K(6, u;)| with \K{b,w)\u{b,w), u{b,w) £ C, then: 

2' = J2 n \^ib, w)\u{b, w) ^±det{K') 

iTi {bw)£m 

where K' : C"^ — > C""' is given by (K'/)(w) — J2be~ ^{b,w)u(b,w)f{b). The undetermined sign is the 
same as before, and thus: 

E( n u{b,w)) = ^ = ^^=dct{K'K-')=dct{K'K-') 

where {K'f){w) = ^bes^ K{b,w)u{b,w)f{b) 

Let us now discuss the height function, following [531 HH] ■ The height function is defined on the "dual" 
"E.^ , defined here as the subgraph of M''' whose vertices correspond to faces of AI that are adjacent to vertices 
of S; we can take the vertices of to be the midpoints of edges of A. Given a perfect matching m of the 
bipartite graph S, one can define a closed 1-form on its dual (ie. an antisymmetric function on oriented 
edges of S^') as follows. Consider the 1-form w G C^(S^) defined by: 



w((6w)^) 



1 if (bw) matched, oriented from black to white 
otherwise 



Then dw e C^(^^) ~ C°(S) is 1 on black vertices, —1 on white vertices, where {du!){f) is the sum of uj on 
the edges bounding the face / of S''', oriented counterclockwise. Given the embedding of S in the plane, one 
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can construct a fixed 1-form ojq £ C^{^^) as follows (here {bw) is an edge of S; the black neighbours of w 
are x, y, x' , y' in this order): 

a;o((5u;)t) = — arg ^^^^ = K{w,b)K-\b,w) "^^ p{w,b) 
2tt y — X 

by Theorem [T] Then duio = duj for any to defined from a perfect matching as above. This is clear from the 
local geometry. 

Since diuiQ — w) = 0, we can write 

ujq — uj = dh 

for some function h on "E) , which is uniquely defined up to an additive constant. In the case of graphs 
in multiply connected domains and on a torus, h is additively multivalued (when tracing h along a non- 
contractible cycle). 

We now seek an interpretation of perturbed operators in terms of height functions. Define 

Ka{w,b) = K{w,b) exp{2i'^ / a) 



Ka{w,b) = K{w, b) exp(23fi ( a) 

J w 



where a is a smooth (0, l)-form: a = a{z)dz. 
As earlier we get 



E|exp(2 J2 ^ [ a)] =det(K„K-i; 
E|exp(2i J2 ^ [ a)| =det(if„K-i 



In the simply-connected case (ie when h is single valued), one can write: 

iMem = ojiiffj) = h{f) - hif) + Miff)) 

where (//') = (bw)^ . Then 



bw)em /est 



{bu 

where df is the boundary of the face / (taken counterclockwise) and P{a) is the R-linear form: 



.6 

P(a) = p{w,b)^ / a 



(bw)eEB 



Let us extend h (initially defined on ^^) to a piecewise constant function, ie constant in each face of Then 
by Stokes' formula 

J2 hif)i^ J a)^^J^hda 



/e3t 
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(here S designates the domain of C covered by the graph; da — da as a is a (0, l)-forni). Finally we obtain 
the following expression for the Laplace and Fourier functional of the height field: 



E ^exp(23?^ hda)^ = dct(K„K-i) cxp(-2F(a)) 
E ^exp(2i3^ hda)^ = det(KaK"^) exp(2iP(ia)) 



The characteristic functional is bounded and is more practical is the absence of boundary. The Laplace 
functional preserves the real structure (Kq, : M"^ — >■ M""'), which is useful for domains with boundary. 

4 Planar graphs 

In order to avoid complications related to non-trivial homology or boundaries, we now discuss the case of 
the plane. 

Let A = A5 be a lozenge tiling of the complex plane C with edge length 6, with S going to zero along 
some sequence; F, is the pair of associated isoradial graphs. We assume that condition (4) is satisfied for 
all S (for a fixed 770). 

From [211 [311 E], we know that there is a Gibbs measure on perfect matchings of S such that, for any 
finite subset {{bi, Wi), 1 < i < n} of edges of S, 

i,wi)em,...,{bn,Wn)em)=\Y\K{wi,bi)] det {K~^{bi,Wj)) 

(manifestly, these local statistics completely specify the measure). In the case where A is the local limit of 
a sequence of biperiodic lattices, the measure can be realised as the limit of uniform dimer covers of (finite 
volume) toroidal graphs. 

Thus we have a probability measure on matchings of S = S^-; E = E5 is the expectation under this 
measure, h the height function, seen as a function constant on faces of S and well-defined modulo a global 
additive constant. 

Consider g G C^(C) a function with compact support, continuous first derivatives and second derivatives 
in Lj^^; set a = dg; we consider again the perturbation 

Kaiwjb) — K{w,b) exp{2i'^ / a) 



Since g has compact support, i^T^K ^ is a finite rank perturbation of the identity, and consequently det{KaK ^) 
is well-defined as a Fredholm determinant. 

Lemma 2. Let h be the height function (constant on faces); then 

det(i4:„K"^) =E(cxp(2i3 [ /laa^)) exp(2i V p(w, 6)5 / a) 

= E(exp(-i5R / h{/Sg)dA))exp{~2iP{ia)) 
Jc 

Proof. We cannot apply directly the finite volume identity we observed earlier. Let us consider v a finitely 
supported function on edges of S. Then 

Yl (1 + w(w,6)l(hu,)em) = X! n l(bt")emi'(w,6) 
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(with finitely many nonzero summands on the right hand side and by convention J^g = 1) and thus 

\{bw)eEs J S={{biWi),l<i<n}CE^ \i=l J 

On the other hand, if K' : C"^ — >■ C"""' is given by its matrix elements K'{w,b) = {1 + v{w , b))K{w , b)) , the 
Fredholm expansion (eg [41]) reads: 

det(Id+(if' -X)K-i) = y det {{{K' - K)K^^){w„Wj)) 

The Cauchy-Binet formula gives 

det m' -K)K-')iw,,w,))= .det UK' ~ K)iwM) det {K-\bk,w,)) 

l<2,7<n ^ — ^ l<t.k<n l<fc,i<n 

~ ~ {bk,l<k<n}cEB ^ ' ~ ~ ~ 

and besides 

71 

det {{K' - K){w^,bk)) ^ V sgn{a)Y[ K{w^,b„(^,■))v{wi,b„(^,^) 

~ ~ crGe„ 2=1 

This completes the identification (reordering rows of deti<ij<„(K~^(6i, Wj)) absorbs sgn((T)): 

El II (l + ^;(^z;,6)l(fc^)e„^) =det(i^'K-i) 
\{bw)eEs / 

Specialising to K' — concludes (the rest of the argument being as in the finite volume case). □ 

The next step is to construct and estimate (in particular near the diagonal) a kernel Sq inverting Ka- 

Let us first consider K^, as a finite difference operator. For (j) a function on C, we define restriction 
operators Rb,Rb by: (i?ij0)(6) = (i?B</')(&) = 0(6) for 6 € F and (i?B0)(6) = -(RB(f){b) = i(f){b) for 6 e F^. 
Let us also denote {RwP){w) = 6(w)e-*''('") for (3 = b{z)dz and {RwP) = b{w)e"'^'"'> for (3 = b{z)dz. These 
restriction operators are C-linear. 

Let w e Mw be an edge node, with black neighbours x, y, a;', y' in ccwise order labelled in such a way 
that [xx') is an oriented edge of F. We have 

K(i?B0) = 2^i^Rw{^4>) + 0((5V (4-2) 
= 2^i^Rw{ddp) + 0((53w^,,(,5)) (4.3) 

where w designates the modulus of continuity, 0' and 0" the gradient and Hessian of 0, and /i^ is seen as a 
diagonal operator on C*^"'. 

Observe that 

K^{Rb^){,w) ^\i{RBm{w) 
where 0(z) — (f){z) exp(2i5 a). We deduce (here a — g^) 

K^{Rbc^) = 2^i^Rw{d^ + 0«) + S^Oiu:^,{S) + SW\UHoo + ll0lU^a('5) + |10|loo'5|la|lL) (4-4) 
= 2^i^Rw{d^ + 0a) + O(53||0||c^(l + ||a||cO') (4-5) 
= 2/i^i?^(90 + 0a) + O(^*||0||c3(l + Mc^f) (4.6) 
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Similarly, we obtain 



Fix Wo G Mw and set A = a(u'o). Then: 

K^{w, h) = e^^'^'^^^^Kiw, &)e-2'»(^^) + 0{5^uja{\h - u^oD) (4.7) 

The finite difference lead us to think of as a finite difference version of (simultaneously) d -\- a and 
the adjoint operator d — a. To construct a kernel Sq inverting Ka, we are going to construct an approximate 
kernel Sq using at large scale inverting kernels for these continuous Cauchy-Riemann operators, at small 
scale {\b — u;| ^ 1) discrete holomorphic functions; and finally control the error between Sa and Sq. 

We have a = dg. Thus d + a = e^^{d)e^ , and consequently 

e9(«')-9(z) 

Sa{z,w) 



7r(z — w) 



is a kernel inverting {d+a) (on the right): {d+a)Sai-, w) — Swdz as distributions. It is uniquely characterized 
by Sa{z,w) — as z -> oo. 

Set 

S„(6, w) = e2'^-(^(''-'"»(K-i(6, w) + Rsifi) + Rb{}^')) 

for 1 6 — I < r/ and 

Sa(6,^) = ^ (i?B(e*''("')5a(6, «;)) + i?B(e-^''(^")5-„(6, zi;))) 

for |6 — iy| > 77, where ry, /i, /i' are parameters to be specified. 
We have 

= ^ . - -mw) + g.{w)) + 0{b -w)= ^ - -{dM9){w) + 0(e2"ff"='c.g,(|6 - 
7r(o— luj TT ■n[o — w) TT 

We can estimate KoSai-^w). If w' is such that all its black neighbours are in 0(10,1]), given that 
K^^{b,w) = 0{\b-w\-^), we get from (|4j|: 

iK^U,w))iw') - SUw') + O ( s^ -A\u''-w\ +S) 

\ \w' — w\ + 



If all black neighbours of w' are outside D{w,ri), we have as in (4.4) 

S , S\\g'\\ 



iK^U;w))iw') ^ S-'e'M-O f + + — 

and for w' outside of D{w,ri) U supp{g), 

iKj^i.,w))iw') = ( ) 



'l|2 

00 



ii^P \w' — w\ 



We have seen that 



K-\b,w) = -f-— + 4== +0 



2 \TT{b-w)J 2 \7r(6-w)y Vl^-^i'l 
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Set 

/i=ie-(-)(--(9,5R5)H) 

I TT 

Then if dist(6, w) is of order t], the difference between the short and long distance definitions of is 
of order 0[ijjgi{r\)e^^3^°° + S/rf'). Consequently, if w' has neighbours both inside and outside of D{w,ri), 

Assume now that g is in the Sobolev space W^'P with p > 2 and has support in -6(0, r). In what follows 
constants may depend on p,r. By Morrey's inequality (eg g' is e-H61der with Holder norm less than 
c||5"||p for e 1 - 2/p. It follows that ujg,{s) < c||.g"||pS^ and \\g\\ci < c\\g"\\p. Set C = ||g"||p. 

Thus, the norm of Ka^^i-i w) — 6^ (wrt counting measure) is less than: 
hence less than 

ce^^(?7^+i + ^ + + d\logr,\ + - + S^) 
T] rj 

Here we simply collect errors stemming from: replacing Ka with an operator conjugate to K at short distance; 
gluing the short and long distance approximation; and using the limiting continuous kernel at long distance. 

Setting now rj — with /3 = '^^^ ^^'^ T = Id —KaSa, we get that 

\\T\\l^ < ce=^((5("+i)/(^+2) v^^) 

where is the — > operator norm, and consequently (Id— T) : ^ is invertible for 6 small 

enough, and we may set: 

oo 

SO that KaSa — Id. 

We now want to estimate ||SQT||/,i_j.ioo . We simply expand 

(SqT)(6, w) = y^^Sa{b,w')T{w',w) 

w' 



and as before we split T in a short range and long range part and notice that Sa{b, w) — 0(6^*^/ dist(6, w)) 
to obtain: 




< ce'^ (77" + (77" + S^-^P)\\og{S)\ + {S'\ log(5)| + S^-^ + 5'-^ + 5^-^^) + J^) = 0(e^^5"') 

The various terms correspond to the possible relative positions of 6, w', w. Notice however that the estimate 
is simpler when b ^ w, which is the most useful case. 



17 



Since - S„ (S„T)(Id -T)"!, ||(Id-r)-i|lii = 0(1), \\SaT\\Li^L-- = 0(e='^5^'), we conclude that 

S„(6,u;) - Sa(&, w) = 0(e'=llf"ll'',5"') 
uniformly in b, w (for p, r fixed; e' > depends on p > 2). In particular for b ^ w we have 

S„(5,w) -e2''-*^(''-"')K-i(6,w) = i^{e"'r^) + Oie''\\s"\\,s<'') 
where = —^{dz^g){w) and ly = u{b) + v{w). 
Lemma 3. The following estimate holds 

\ogdei{K^K-^)^2iY,p{w,h)^ j j |V5R5|2 + 0(5^') 

Proof. Observe that Ka\^^ and KSq, are finite rank (hence trace class) perturbations of the identity, and 
(ifaK"^)(KSa) = Id (since K"^ is also a left inverse of K). If a depends smoothly on a parameter t (say 
a{t) — to), we have the variational formula (eg [20], IV. 1) 

^logdet(i^„S) = Tr((^7^,)S„) 

as long as K^K^^ is invertible, which is at least true for small enough t. From 

Ka{w,b)^K{w,b)exp{2i^ f a)2i^ 

J W 

S^ib, w) = exp(-2i3 / a)\<r\b, w) + i5(e*(''W+''WV„) + 0{5^ ) 



(for b ^ w), we obtain (if x,y,x',y' denote the black neighbours of w in ccwise order) 

K^Sa{w,w) = 2iJ2piw^b)'^ f a-\y' -y\.\x' -x\ (3(Ae-*''("'))5(e''^("'V„) + !3(Ae-*''('")(-i))^J(e'''^"'^i^«)) +0{S^+'') 

b^w Jw ^ 

= 2i^p{w,b)^ / d + 2/io(u')5R(Ar„) + 0((5^+^ ) 

I, Jw 



b^w 



and 



d 
di 



logdet(KaK"^) = 2i ^p(w,6)3 / d + 5 / a A {r^dz) + OiS"') 

b^w 

(since dz A dz — —2idA). Moreover, with a = dg, 

5 J aA{r^dz) = -^^ J dgAd{^g) = ^^ J gOdi^ig) = J g^gdA^-^jJ |V3?5|' 

In particular, for 5 small enough dei{Ka\^^^) does not vanish along the interpolation path, thus (eg [4T1 
chapter 3) i^T^K"^ stays invertible and the variational formula is legitimate. □ 

From the previous two lemmas we immediately conclude: 

Corollary 4. If g G W"^'^ H C^, p > 2, then J converges in distribution to a centered normal variable 
with variance ^ J iVgp, as S '\ 0. 

Thus we can integrate the discrete height field h against a test function in L^, p > 2 (the compact support 
assumption is for technical convenience). An optimal statement (given the free field limit) would involve a 
test function in H^^; however h (as defined here) is not in H^, due to jump discontinuities on edges of S. 
It is unclear whether one can define a better interpolation of the discrete function for which one could relax 
substantially the condition g G W^'^ (to g E say). 
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5 Toroidal graphs 



As before, we start froni a rhombus tihng A = A5 of the plane with edge length i5 ^ 1, from which we 
construct r, rt, = At, M. Additionally, we assume that these structures are biperiodic, with periods 
1 and T = Tg (9r > 0). Denote by T the lattice Z + tZ; by quotienting, we obtain a rhombus tiling of 
S = C/T. We are interested in perfect matchings of S = M/T. Notice that Euler's formula applied to F 
shows that = \'5w\- 

The height "function" h is associated to a perfect matching m on S by a local rule. Since the torus is non 
contractible, the height function is now additively multivalued. The current J = dh is a, well-defined closed 
1-form on S (for some extension of the discrete height function on to a piecewise continuous function on 
S; then J is a distributional 1-form). The Hodge decomposition of J reads: 

J = ujh + dho 

where Wh is a closed harmonic form (hence, in flat metric, coh = adx + bdy for some a, 6 G K) and /iq is a 
(single- valued) function on E, well-defined up to an additive constant. The "topological" component wq is 
uniquely specified by its periods, which equal those of the current: uJh = J J, where 7 is a basic cycle on 
S. 

We are interested in the asymptotic distribution of J, or equivalently the asymptotic joint distribution 
of {uh, ho). For clarity we will be discussing first the topological component of the current. 



5.1 Flat line bundles 

For A e C, denote: 



Kx{w, b) = K{w, b) exp(2i9 / Xdz) = K{w, b) exp(2i3A(6 - w)) 

J w 

so that K\ defines an operator C^® — )• and by quotienting an operator C""' . We are concerned 

with the inverse of that last operator. Let x • ""^ ~^ U = {2; € C : |2;| = 1} be a unitary character; consider 
the finite dimensional space 



:^-)^ = {/ e C^- :yzeMB,0JG T,f{z + uj) = x{'^)f{z)} C C^- 



and 



'^"')x defined similarly; note that 



)ld^ 



(here Id denotes the trivial character). Plainly, 
(T,/)(2)=/(2)e2*^^^ 



as Kx commutes with the action of T by translation, it defines an operator Kx : (C^^^),^ — > (C^*^);^. Set : 



If xa(w) = e^'®"^" is the character associated to A , maps (C*^ )id to {C'^ )-)^^^, and we have the diagram 



K 



XX 



(C^w) 



XX 



(C^«) 



Id 



Kx 



Y 

(C^'^-)id 



so that we may focus on K : (C*^^);^ — >■ (C''^'^);^. An element in the kernel restricts to a bounded harmonic 
function on F and F^'; then (Liouville) these restrictions are constant, hence the kernel is trivial ifi^ X 7^ Id. 
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As |Sb| = |Sm/|, K : (C^^«)x ^ (C^"')^ is invertible iS x W- Let denote the inverse operator (x ^ Id). 
We now want to relate S^^ to kernels inverting (continuous) Cauchy-Riemann operators. 

Denote by the holomorphic, unitary line bundle over E obtained by twisting the trivial line bundle 
^id by the unitary character x ■ ~^ U. Sections of may be seen as multiplicatively multi- valued 

functions on S, with the transformation rule prescribed by x- + •2^) = xl"^)*!^)- 

We may consider the operator 

canonically associated to the complex structure on L^. It is well known that this is a zero index operator 
which is invertible iff x 7^ Id. Moreover its inverse has an explicit expression in terms of ^-functions, which 
we now recall (see eg [IS])- 

Consider the 6 function with characteristics: 



1? 



2e 
2e' 



(z) = ^exp2i7r ( ^T(n + e)^ + (n + e){z + e'U = exp(2i7r(^£^ + ez + £e'))-d 



(z + e' + £t) 



which transforms as: 



2e 
2e' 

2e 
2e' 



(z + 1) = cxp(2«7r£)i? 



2e 
2e' 



{A 



(z + r) = exp(-2i7r(z + e' + -))i9 



2£ 
2£' 



(^) 



2£ 
2£' 



(see eg [TS]). For concision, let us denote i} = i} 
are as in [5]). Then consider the meromorphic function 



and 6* = i9 



(these last notations 



Then T(z + 1) = -e^'^^T^z), T{z + r) 
residue 1 (if (2£,2£') ^ (1, 1)). Thus: 



^^^^ d{Q)e(z) 

_g-2i7re xi^z). Moreover T{z) has a simple pole at z = with 



S^{z,w) 



1 d{z~w)e'{Q) 
^ ' i?(O)0(z- w) 



is the kernel inverting 9^, where x is the character given by x(l) = — e , x{t) 
may take 3A — 7r(£ +5), 3(Af) = — 7r(£' + \), so that 



, X 7^ Id. We 



A 



r + 1 



£ + £T ■ 



Since is odd, we have the asymptotic expansion near the diagonal: 

1 



S'x(z,w) 



7r(z — w) 



+ r^{w) + 0{\z-w\/d{Q)) 



as z — > w, where 



Ttr^iw) = lim {■kS-^{z — w) 



^)4(0)=2^.£+|(£'+£.) 



(which actually depends only on x)- If we set = 2 — x(l) — x(''')i which is comparable to (e — 1/2)^ 

(£' - 1/2)2 foj. g [Q^ J^^yg ^(-Q) ^ 0{^). 
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If £, e' depend smoothly on a parameter u, we have: 



d_ 
du 



log I? 



2e 
2e' 



(0) ] = (iTrre^ + iiiree') + log 6I3 (e' + er) 
/ du au 

= 2i'Ke{iT + e') + 2i7ree' + (e' + iT)-^{e' + er) 

= (A3T)rx(0) +2i7ree' 

the last term being pure imaginary. 

Following Ray and Singer f|39j). we consider the analytic torsion defined by 

T^ix) = exp(-ic;(0)) = det^{idj*d^y/' 

where Cx is the C function defined from the Laplacian-type operator {d^)*d^. Then an explicit diagonalisation 
of {d^)*d^ (in flat metric) and Kronecker's second limit formula implies that ([32]) 



-\e{li-Ty)/7^{T)\ 



where x(mr + n) = exp(2i7r(r7i/i + ni')) and 77 is the Dedekind rj function (eg [5]). Here /i = |— £', v = e—}^, 
so that 



rE(x) = e--(^^)'/'^n^3(^)M^)l = l^ 



2e 

2£' 



(0)/r/(r)| 



with z = e' + er. Note that 



2e 
2e' 



(0) 



-7r(9z)^/QT| 



j(z)|, where z = e' + er. 



We admit for now the following near diagonal estimate for S^: 

S^{h,w) - K~^(6, w) = i9(e*'^rx(0)) + 0{5^) 

for |6 — w| = 0{5) (in particular if 6 ^ w in A/), for some positive constant e. We also admit that the 
estimate is uniform in x foi' X in a compact set not containing the identity. 

Let us consider now A as a differentiable function of a parameter u and analyse the variation of det K\ , 
Kx ■■ C"^'' . Assume that x\ 7^ W, so that Kx is invertible. Then: 



du 



det{Kx)^TTikxK^') 



where for b ^ w, 



kx{w, h) = 2i5(A(& - w;))e2*^(^(''-"'»K(6, w) 



= K-\b, w) + i^e'^'r^iO)) + 0{S/c^) 

where x = X\- If ^ in F (resp. F^), w the white vertex corresponding to {bb'), we get K(w, &)K"^(6, w) = 
K(u', 6')li^^(^': ^) f^rom Theorem[l] Thus the contribution of K~^(&, w) in the trace cancels exactly; besides, 



2i'^{\{b - w)))K{b,w){i'^ie"'r^iO))) = -^i<:^iw)'^iXe-"')'^ie'''r^{0)) 
(where v — 1/(10) + t^(fo)) and we are left with: 

Ti{kxK^^) = 2 ^ M<>(w)3?(Ar^(0)) + 0{S/c^) = 2.Area{C/T)di{Xr^{0)) + 0{S') 



23?^ log 
du 



2e 
2e' 



(0) 



0((5^) = 2— logTs(x)+0(5^) 
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since around alternatively 1 and i; and Area{C/T) ~ Sr. We conclude (admitting the near 

diagonal estimate) that if O Xi noii trivial (i = 1,2), then 



as the mesh (5 \ 0. 

5.2 Bosonisation identity 

Let (j)nm be the harmonic differential on S with half-integer periods n — (f'nm and m = Jg 4>nm- Let 

S{(t)) = 2ti / \V(t>rdA + 3((m - fn)z) + Amnm 

= — |m — fnp + 4i7r(m£ + ne') + Amnm 

where z = e' + st = x + iy, and 

Then a Poisson summation argument ([31, Section 4.C) shows that 

Z,„,, = (23r)5e-2^(^^)'/'-n^3Wl' 

The determinants det(iirA) count dimer configurations with some unitary weight, which depends only on 
the periods of the current. This is originally due to Kasteleyn ([23]); see [5] for a recent (and exhaustive) 
treatment. 

We begin with det K (which is as the kernel of K contains constant functions) . Then (for a proper 
choice of ordering of vertices), each term in the determinant expansion corresponds to a matching of S 
(with associated current J) counted with a positive sign if ( J, Jg J) — (0, 0) mod 2, and negative sign 
otherwise. Let us denote this sign by Q(m) = Q{J). 

In the expansion of det-fCA, each matching is counted with an additional phase: 

cb 



a 



ibw)e 



where a — Xdz. As before, we can integrate by parts over a fundamental domain C bounded by cycles A, B 
drawn on F to obtain: 

E E (M/')-M/)-^o((//'))) r« 

dh a — / dh a 
B J A J A Jb 

(Note that in this case, E'^o((//')) K" " = 0). Thus: 

Z-'^detiKx) =E(^Q{J)exp{2iQ{-Xf J J + J J)) 
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where Z is the partition function of the dimcr model on S. Let us set (n, m) = ^(/^ J, /g J), the half-periods 
of the current, and 

A = A £,£' = ^ U + 



as before. Then: 

Let us remark that 
so that 



2i3(-Af / J + A / J) = 4i7r(me + ne') + 2i7r(r7i + n) 
J A Jb 



Q{J) = -e 

det(ii:A) = E /'e4i'^('"^+"'^')+4»7rmn 



which is 1-periodic in e,e' . 

Let us admit for the time being that 



deti^As f Ty^{X2)\'^ _ Zinst{z2) 



deti^Ai \Tt.{xi) J Zinst{zi) 
for Xi, X2 7^ Id (as we have seen, this follows from near diagonal estimates). 

This is enough to conclude that (m, n) converges in distribution. Indeed, if we write K^^g,/ = Kx[e,E')j 
have: 

(detKe,£' + detK^+i,,, + detX^ i - det i i ) = 2E [e^^^(rne+ne')\^ 

by virtue of 

for m,n E ^Z. If (e, e') = (i, i), det i4r(£ g/) = det K = (as constant functions are in the kernel of K) and 
Zinsti^' + £t) — 0. Consequently, 

^-^(detXo.o + deti^i + detXo i) 2E(1) = 2 

or my- 

Z = - (det Kofi + det -fC i o + det Kq i ) 

^ 2 ' '2 

It follows that, for any (e,e') G (distinguishing the case (e,e') G (|2)^, where we use det(K) = 0), we 
have 



E 



^g4i7r(m£+ne')^ ^ ^g4i7r(me+ne')^ 



where Eq is relative to the probability measure on with weights proportional to e i'l™ "^"1 . Con- 

vergence of the characteristic function then implies convergence in distribution of the half-periods (n,m). 

5.3 General Cauchy-Riemann operators 

Let a — a(z)dz be an T-periodic (0, l)-form (equivalently, a (0, l)-form on E). We can define a perturbed 
operator 



so that Ka defines an operator C*^^ C^^"' and by quotienting an operator C"^ C"** . We are concerned 
with the inverse of that last operator. We first consider the limiting continuous Cauchy-Riemann operators. 
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The Dolbcault decomposition of a reads: 

a = Xadz + dg 

where g is a function on E, Aq e C constant. Thus d + a — e^^{d + Xodz)e^ . In the case x — XAq 7^ Id, this 
is invertible, with inverting kernel given by 



From the expansion S^{b,w)e'^'^^^°^'' '"^ = l/7r(6 — w) + + 0{\z — w\/^/c^), we obtain the asymptotic 
expansion as & — ^ w, w S S fixed, A = a(w): 

\7r(6 — w) / 

^ +rx- Ug.iw) + 9.{w)) + 0{\b - w\.\\g\\c2/^c^) 



Tl{b ~ w) TT 



Let us turn now to the discrete operator Ka, seen as a finite difference operator. For a function (j) 
on S, we define as before restriction operators Rb,Rb by: {RB(t>){b) = {RB4>){b) = <l>{b) for 6 e F and 
{RB<j>){b) = —{RB(f)){b) — i4>{b) for b G F'''. Let w S be an edge node, with black neighbours x,y,x',y' 
in ccwise order labelled in such a way that {xx') is an oriented edge of F. We denote {Rw(3){w) — 6(w)e^*'^*^'^'-' 
for P = b{z)dz; then (see ( |4.4[ )) 

K^{RB<t>) = 2tJL^Rw{d<t> + ^a) + 0{5^\\<t>\\c^{l + HalbO^) 

Similarly, with {RwP) = fe(w)e*''("') for /3 = 5(z)dz, we have: 

= 2^l^Rw{^^ - cj>a) + 0{6^U\\c-4l + ||a||cO') 

Fix Wo G Mvi^ and set A = a(u'o). Then: 

&) = e2*'^(^'")K(u;, 6)e-2^^(^^) + 0((52 dist(6, u;o)l|a|lcO 



Set 

S^{b,w) = e2'»(^(''-'^»(K-i(6, w) + + RB{^i')) 

for 1 6 — w I < 77 and 

for |6 — > r/, where 77 is a mesoscopic scale and 

Z TT 

Reasoning as in the planar case (it is a bit simpler here as E is compact), we may use this approximate 
inverse to show that is invertible for 5 small enough and its inverse Sq satisfies the following near diagonal 
estimate: 

S„(6, w) - e2*'-*^(^-"')K-^(6,u;) = i'^{e"'r^) + 0((5'^) 

where ^ — ^{dzdi.g){w), e is a positive constant, and the error term is uniform in (Ao,.g) for H^Hc^ 
bounded and x = XAq ^ compact subset of the group of characters not containing the identity (r^ blows 
up as X ^ Id). 
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Similarly to the planar case, this implies: 



detK. 
detK. 



OCl 



exp(2z(P(ia2)-P(iai))) -^exp(^-^ (^JjV^g^l^dA- JjV^g,\'dAj^ (^^Mgj 



where a, = Xidz + dgt, Xi is the associated character, and convergence is uniform for H^iHc^ bounded 
and Xi away from the trivial character. Let us point out at this stage that the additive decomposition 
of Ta reflecting the Dolbeault decomposition of a leads to the multiplicative decomposition of the relative 
determinant det{Ka^) / det{Kai) (in the fine mesh limit), which in turns yields the independence of the 
instanton and scalar components of the limiting field. 

To complete the identification, we notice that 

det(Aa) exp(2iP(ia)) = ZE ^exp(4i7r(TOe + ne') + Ammn) exp(-i3? j h{Ag)dA)^ 

and 

Z = - (det Ko,o + det Kiq + det Kqi) 

which implies 

E ^exp(4i7r(m£: + ne')) exp(— iSR J h{Ag)dA)^ — 5^ Eq (exp(4i7r(me + ne'))) exp j \^^9?dA 

in the generic case ^ (5^)^- In the case g = 0, we already noticed that this ensures convergence in 

distribution of the half-periods (n, m). Then for (no, mo) S (5^)^, we write 



^[0,1]2 



and 



E (l{(n,m)={no,mo)} exp(-iJi ^ /l(Aff)rfA)^ = ^-Ai^(moe+n,e')^ (^^U.{me+n.') exp(-iSR^ K^)dA)^ dede' 

Since P((n, m) = (no, too)) — > Po((n,TO) = (no, too)) > 0, we deduce by dominated convergence: 

E ^exp(-iSR j /i(^)rfA)| (n, to) = (no,TOo)^ — )• exp ^— — j \V^g^dA^ 
for any (no. Too) e (|Z)^- 



^2' 

Prom this we deduce: 



Proposition 5 (Finite dimensional marginals). Let ai,...,ak he l-forms on S. As 8 \i 0, the joint 
distribution of 



Js A *«!, . . . , / Js A *ak 



converges to the joint distribution of 



J J A *«! , . . . , y J A *a/c^ 



where J is the current of the compactified free field on S with compactification radius 1 and action S{J) = 
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Proof. Wc have the Hodge decomposition a.j = {ajdx + bjdy) + dcjj. gj E C^. We know that the marginal 
distribution of the current periods converges, and that the conditional characteristic function 



E(exp(j Xj J J A*dgj)\{n,m)) 



converges pointwise. This is enough to ensure convergence of the joint distribution. □ 

To obtain a functional CLT from this finite dimensional CLT, we need a tightness estimate. 

Lemma 6. For any e > 0, the family of probability measures induced on H~'^~^{fl^{T,)) by the current Js 
is tight. 

Proof. Let us remark that for a r.v. X, 

E(X2) < E(0o,o^') +E(go,iX2) +E(Qi,oX2) 

where Q^.e' = e*''^(™^+"^')+^''^™" (a random sign), since 1 < Qo,o + Qo,i + Qi,o < 3. 
Set a(e,£') {t) = A(e, e')dz + tdg and X = / /lAgirf^. We get 

E((3? j hoAgdAf) < -Z"! — ^^^^(det if„o,o(t) + det + det K«i,o(t)) 

We have: 

detJr„^^,(t) = ZE(g,,e'e-'*^)exp(-2ztP(ia5)) 
and thus ^ 

E{Q,,,,X^) = -2Z-'^ (deti^„ eM'^itPiidg))) 



We now assume that {s,s') G {(0,0), (0, 1), (1,0)}, so that in particular det-K'^,^ is invertible at t = 0, 
and 

I logE(ge,.'e'*^) = Ti{K-^K„) + 2iP{idg) 

^ logE(Q,,,,e^*^) = Tr{K-'K^) - Tv{{K-'K^f) 

where a = ol^.s' {t) for brevity. 

From the short distance asymptotics for K~^, we get 

TriK-^K^) + 2iP{idg) = 0{S'"\\g\\ci) 
at t = for some e" > 0. We also have 

Tr{K-'K„)^0{\\g\\ci) 

We are left with estimating Tt{{K-'^ Ka)"^) (at t ^ 0). Since K-'^{b,w) = K~^{b,w) + 0(1) and Ka = 
0{S^{1 + llslIcO)' by isolating the leading singularity we may write 

Tr{{K-^K^f)= J2 {KaK-'){w,w'){k^K-'){w',w) + 0{{l + \\g\\cr)^ 

\w — w'\<€q 
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where eo < min(l, 9r)/10, say. Fix u>; replacing Ka with where 13 = tgg{w)dz (a constant (0, l)-form) 
induces an error of order 0(||(7||ci,e (1 + ||g||ci) (here C^'^ designates functions with e-Holder derivative, 
e > arbitrarily small but fixed). Thus we simply need to estimate 

w':\w—w'\<eQ 

which may be thought of as a discrete version of a principal value integral of type p. v. J J ^^^_i™,-j*2 ■ A 

probabilistic interpretation of this quantity goes as follows: a product K~\b,w')K-\b' ,w) (for b ^ w, b' ^ 
w') is, up to multiphcative local factors, the covariance Cov(l(f,^)£^, 'i-{h'w')em) under the appropriate Gibbs 
measure on tilings of the full plane. By linearity {w is fixed and we sum over w'), we are left with estimating 
Cov(l 

(bw)£xm^)i where ^ is a linear function of the heights in B{w,eQ). Since /3 is constant, dfi — and it 
follows that £ depends only on the heights on dB{w,eQ). Since Cov(l((,^)£m) l(6'«;')em) = 0{l/\w' — w\^), 
we conclude: 

W' •.\w — w'\<€q 

We finally get the estimate 

E((5R j ho/^dAf) = 0((1 + ||5||c0(l + llsllc^.O) 

which is uniform in 6 for 5 small enough (the limit as i5 \ is of order HffH^i). 
Consider an eigenbasis for the Laplacian on S = C/T: set 

Quiz) = exp(i3fj(zu)) 

where u e T = {t; e C : Vz G T,^{zv) e 2ttZ}. Note that WguWc" = 0(1 + |m|'=) and ||5„||ci.- = 0(1 + |m|^+^). 
We may define 

\\J2^-9u\\is = Y^K\'{i + \u\y 

If we choose (which is given modulo an additive constant) so that J^, hodA = 0, we may write ho = 
E„eT\{o} "-udu, where 

E((a„)2) < c|w|^-2 

for M e T \ {0}, and c is uniform in 6, u. Consequently, the Chebychev inequality yields: 

P(V«eT\{0},|a„|<CH-^)>l-^ J2 H^-'-'^ 

«et\{o} 

where the sum converges if 27 + e < 0. On this event, 

« uet\{o} 

which converges if 2s — 27 < —2. By taking 7 e (— e, — e/2), and observing that H^^{T,) is compactly 
embedded in H^''{T,) for si > S2, we conclude that the probability measures induced by ho on H~^~^{T,) 
arc tight for 6 small enough. 

The current can be decomposed as Jg = dho + (jJh- Since ojh takes values in a two-dimensional lattice 
of H^{Q^{T,)) and converges in distribution, it induces a tight family of probability measures (s arbitrary). 
The lemma follows. 

□ 
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We may now state a functional limit theorem for the current. 



Theorem 7. For any e > 0, the probability measures induced on _ff^^^^(f2^(S)) by the current Jg converge 
as 6 \i to the distribution of the current of the compactified free field on S with compactiflcation radius 1 
and action functional S{J) = ^ J J A * J. 

Proof. We have obtained tightness of the measures and convergence of the characteristic functional. The 
dual of may be identified with H'^'^^ {il^ (Y,)) , which consists of 1-forms with coefficients. 

Indeed, 77^+^ may be embedded in M/2,p+e j^gg jjj^ 7.58), the elements of which have continuous versions with 
e"-H61dcr derivatives (Morrey). For Banach space- valued variables, tightness and pointwise convergence of 
the characteristic functional ensures weak convergence (eg [33j, 0.2.1). □ 

6 Surgery 

In the most general set-up, we shall be interested with correlations of electric and magnetic operators for the 
height function in the plane, on a torus or in a finite planar domains. This corresponds to a certain Cauchy- 
Riemann problem, with singularities at insertions. In this section, we show how one can treat separately 
issues stemming from boundary conditions, electric insertions and magnetic insertions and then glue these 
components together. 

We start from a family of rhombi tiling {As)s with edge mesh 6 going to zero along some sequence. Let 
Di (resp. Do) be a simply connected neigbourhood of (resp. oo) in C that intersect in an annular domain 
A. Let 7 be a simple closed loop such that A is a tubular neighbourhood of 7. Wc assume that the distance 
between 7 and dA is large enough compared with 6. 

Let Si = Si((5) (resp. So) be a bipartite graph that agrees with AI in Do (resp. Di); Ks : — >■ C"= , 
s € {i,o}, is a linear operator such that Ki (resp. Kq) agrees with K in Do (resp. Di); in other words, 
"BijKi are obtained by modifying M, K in 1?^ \ Do, and vice versa). We assume that Kg is invertible in that 
for each w G S^, there is a unique function / G vanishing at infinity such that Kgf = (5^, s G {i,o}; 
it is denoted K~^{.,w). This implies in particular that Kgf — and / vanishes at infinity iff / = 0. We 
assume that, for s G {«,o}, there are kernels {z,w) Ss{z,w), {z,w) 1— !• Ss{z,w), and rj is an error rate: 
lim^x^^o vi^) — such that: 



on A"^ \ Aa, uniformly on compact subsets, where Aa = {{x,y) G A"^ : x y}. Note that we do not request 
that Si,{z,w) = Ss{z,w), which is the case when is real. More precisely, we assume that 



is in A^ and holomorphic in z in A (in the presence of boundaries, it will be harmonic, rather than 
holomorphic, in w). The holomorphicity condition is actually superfluous but will always be obvious in 

applications. Correspondingly {z,w) 1— >■ Ts{z,w) = Ss{z,w) ^ ^ is in (z,w) and antiholomorphic in 

z. 

A discrete Cauchy integral formula argument (see below) shows that the convergence assumption is 
equivalent to assuming: 



K-\b, w) = (£'''('") 5. (6, w)) + -RBie-'^^^^S^ib, w)) + 0{r^{6)) 



{z, w) i-> Ts{z, w) = Ss{z, w) 



1 



7r(z — w) 



K-\b,w)^K-\b,w) + -Reie' 



'^^^Ts{b,w)) + -RB{e-'-^'"^f,{b,w)) + 0{v{S)) 




uniformly on compact subsets (including on the diagonal). 
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We now consider the glued data: Eg agrees with in D^, s € {i,o}; Kg : C~s C~s agrees with 
in Ds, s S {i, o}. We wish to estimate (if defined). 

Let us consider 7^ a simple cycle on F (and thus on M) which approximates 7 in the sense that each arc 
of 7 of length £ is at Hausdorff distance < C6 of an arc of 75 of length at most Ci, C > fixed. 

Let / e C~* be such that Kif{w) = for any w e which is (strictly) inside 75. Let f{b) = f{b) if 
b is on or inside 75 and otherwise. Then Kif is supported on white vertices on 75 or adjacent to a vertex 
on 75. Moreover / — K~^{Kif) = (since it is in the kernel of Ki and vanishes at infinity). This yields the 
Cauchy integral formula: 

m = Y,K7\b,w){Kif){w) 

w 

for b on or within 75. Let 7^ be the set of black vertices which are either on 75 (and thus on F) or inside 
of 75 and adjacent to a white vertex on 7^ (and thus on F^). Correspondingly, let 7^ be the set of white 

vertices which are either on 75 or outside of 75 and adjacent to a black vertex on 75. Let K^{w, b) = K(w, b) 
if w G 7*^, b G and Ky{w, b) = otherwise. We can rephrase the Cauchy formula as: 

m= Kr\b,w)iK^f\,B){w) 

for b G on or inside 75 and / such that {Kif){w) = for w inside 75. The (discrete) Cauchy data space 
(on 7, for Ki) is the subspace Cf C C"^ consisting of restrictions to of functions / such that Kif = 
strictly inside 75. Clearly, 

pf : C'^'' — > C^" 



Kr\.,w){K,g){w) 



is a projector onto the Cauchy data space . We want to relate this to limiting continuous Cauchy data 
spaces. Given our data, we can define them in terms of the following operators on functions on 7: 



iw 



Pi{f){zo) = \im ^ (f Si{z,w)f{w)d'h 
Pi{f){zo) = lim f Si{z,w)f{w)dw 

z^zo Zl 

where limits are taken from inside 7. Writing 

^ j( Si{z, w)f{w)dw =^.j> Ti{z, w)f{w)dw + f{z) + j( ^^^^0^dw 

gives that Pi is a bounded operator Lip(7) — >■ C°(7). This defines Cauchy data spaces by Ci = {/ G Lip(7) : 
/ = Pif}, Ci = {f £ Lip(7) : / = Pif} (more classically one would consider the L"^ closure of these; 
Lipschitz functions are sufficient for our purposes). 

Let / be Lipschitz in a neigbourhood of 7. We wish to estimate the discrete "contour integral" 

^ Kr\b,w){K^iRBmw) 

where (J?b/)(6) = f{b) if 6 e F n 7^ and {RBf){b) = if{b) if 6 e Ft n 7^. For this we note that 7^ 

consists of white vertices on 75 (a simple cycle on F) and white vertices on 7] (a cycle on F^), The path 7J 
may be described as follows. Let (60, 61, . . . , 6„ = 60) be the black vertices on 75 (taken counterclockwise). 
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For each i, enumerate (in counterclockwise order) the faces of T which are adjacent to bi and outside of js- 
^Ii'---'^lfe,- Concatenating these hsts, one gets 7] = (bj . . . , 6| 1, . . . , , . . . , fcj^ 1 = feji), a 

cycle on (which may involve some backtracking). Without loss of generality, one may assume that the 
reference orientation (used to define the Kasteleyn orientation) of edges of F on 75 agrees with the direct 
orientation of 75. Then if w on 75 corresponds to the edge {bb') of F, 

{K^iRsmw) =^-\b'- b\im + 0{6'\\f\\Lip) 

and if w on 7] corresponds to the edge {bb') of F^, 

{K^iRnfrnw) = sgniKibo,w))^\b' - 6|/(6) + 0{S^\\f\\up) 

and bo is the black neighbour of w which is on 7''' (if w has two neighbours on 7^, {Kj{RBf)){w) = 
OiS^WfWup)). Then 

^ ie'-'^^^T,{z,w)){K^iRBmw) = l(f Ti{z,w)f{w)dw+l<f T,iz,w)f{w)dw + 0{Smup\\Si\U 

and by Stokes' formula, Ti{z,w)f{w)dw — Ti{z,w)f{w)dw = 0((5||/||iip||Ti||oo) (as the area of the 
annulus between 75,7] is 0(1))- (Estimates are uniform for 2: in a compact subset of A). Similarly, 





'"'^'"^Tiiz, 


w)){K^{RBmw) 










"'^'"^Ti{z, 


w)){K^(R^f)){w) 










'•'^'"'^fiiz, 


w)){K^{RBmw) 








E(- 


''^("')Ti(^, 


w)){K^(R^mw) 









= -i(f fi{z,w)f{w)dw + 0{d\\f\\up\\Ti\\oo) 



To deal with the singular part, we observe that the constant functions RBin), Rb{ij) are discrete holomorphic 
and consequently by replication: 

\k-\b,w){K^{RB{ij))){w) = RBilA-^^h) 

E \<.-\h,w){K^{R^{,j))){w)=R^{y)l^.{b) 

where 7' is the set of vertices on or inside 75. Then 

K-\b,w){K^iRBmw) = RBUmi^^ + Y K-\b,w){K^{RB{f - f{b)))H 



RBimn^ + ^Rb ^^'^^^_j^^^^ d^L^+0{5\\og5\.\\f\\Lip) 



and similarly 



Y ^-\b,w){K,{RBmw) = RBimn-y. - ^Rb U l^^^^l^Mdw) +0{S\log5\.\\f\\Lip) 
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We conclude that if / is Lipschitz around 7 (||/||lj:p its Lipschitz norm in a neighbourhood of 7), then 

\\Pf{RB{f)) - RB{P^f)\\o. = Oif,i6)\\f\\L^p) 

(on 7), where ry((5) = ri{S) V ((5|log(5|). In particular, if / is Lipschitz around 7, then it is in the (con- 
tinuous) Cauchy data space Ci iff ||P/(i?B(/)) — i?s(/)||oo goes to zero as (5 \ (uniform norm on 7-^). 
Correspondingly, with the same assumptions, 

\\p■(R^{f))-R^{p^f)\\o. = oifiismu.p) 

again on 7. 

Having analysed inner Cauchy data spaces, one may repeat the argument for outer Cauchy data spaces. 
Let 7^ be the black vertices which are on 7a or are outside of 75 and adjacent to a white vertex of 75. The 
outer Cauchy data space C C"" consists of restrictions to 7^ of functions / G C"° such that Kof{w) = 
for any w G 2^ strictly outside of 75 and / vanishes at infinity. The continuous outer Cauchy data space 
Co, Co are defined as fixed points of the operators Po,Po- 

PoU){zo) ^ \im ^ (f So{z,w)f{w)dw 
Po{f){zo) = lim (j) So{z,w)f{w)dw 

z-i-zo 2i 

where limits are taken from outside 7. 

Let us address uniqueness for Kg, ie Kgf = 0, / vanishing at infinity implies that f ~ 0- This will follow 
(for small enough d) from the following natural assumption on continuous Cauchy data spaces: C^nCo = {0}, 
Ci n Co = {0}. Indeed, assume by contradiction that for some sequence (5„ \ 0, there is /„ G C"9 " 7^ 
such that Kg" fn — 0, /„ vanishes at infinity (the line of argument here is similar to some arguments in 
[?])• Let 7 be a simple cycle in A that disconnects from 00. We normalise /„ so that ||(/ri)|'y||oo = 1- 
Let 71,72 be two disjoint cycles in A bounding an open annulus A' that contains 7. By replication, we 
see that ||(/n)|7iU72lloo is bounded and consequently the Lipschitz norm of /„ on compact subsets of A' is 
bounded. Up to extracting a subsequence, one may assume that (a suitable interpolation of) (/„) converges 
uniformly on compact subsets of A' to a non vanishing (since it has uniform norm 1 on 7) Lipschitz function 
/. More precisely, we may write /„ = i?_B(ffn) + Rsihn), where ((/„) and (/i„) converge in Lispchitz norm 
on compact subsets of A' (again by replication). Consequently, g ~ lim is in Ci H Cq = {0}, /i = lim /i„ is 
in Ci DCq — {0}, yielding the needed contradiction. 

We may now address the problem of gluing inverting kernels. Assume that Sg : ^ C is a kernel with 
the same regularity conditions as Si, So and such that: for any simple cycle 7 in j4 disconnecting from 
00, if w is outside 7, Sg{.,w) G Ci and Sg{.,w) — So{-,w) e Co ; and if w is inside 7, Sg{.,w) G Co and 
Sg{., w) — Si{., w) G Ci (where Ci, Co are the inner and outer Cauchy data space on 7). Similarly, we assume 
given Sg : ^ C a, kernel compatible in the same way with inner and outer Cauchy data spaces Ci, Co, and 
with the same regularity as Si, So- Given this data, we start with constructing Sg, an approximate inverting 
kernel for Kg : C"s — > C"9 (at least for some w's). 

First let us consider two disjoint simple cycles 71, 72 in A that disconnect from c», with 71 inside 72; 
7* is an approximation of 7i on Ts. For w € within 0{S) of 72, we set: 

Sgib,w) = K-\b,w) + pf^ (^^RBie''''^'"\Sg{.,w) - Soi.,w))) + ^i?^ (e-'^W (5,(., «;) - Soi.,w))))^ 
for b outside 71 
= Pf^ (e-W5,(.,u;)) + ^i!^(e— W5,(.,u;))) 

for b inside 71 
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and for b on 7^ one may use either definition (here P^^ , PJ'' denote the inner and onter discrete Cauchy 
data space projectors described earUer). Symmetrically, for w G within 0{6) of 71, we set: 

S,{b,w) = K7\b,w)+Pf^ (^^RB{e"'^'^HS,{.,w) - Si{.,w))) + bis [e-"'^^\S,{.,w) - «;))))) 
for b inside 72 
= Pj^ (e-W5,(.,u.)) + ii?^(e— 

for b outside 72 

Let us observe that if w is within 0{5) of -ji, KgSg{.,w) — Sy, is supported on white vertices with graph 
distance < 1 to 7f_j, i € {1, 2}; let us denote 7]^ these sets of white vertices, i G {1, 2}. 

In the continuous limit, the assumptions on compatibility of Sg with the Cauchy data spaces specified 
by Si, So translate into the replication identities: 

P;''{Sg{.,w))=Sgi.,W), P^^{Sgi.,w) - Soi.,w)) = Sg{.,w) - So{.,w) 

on 71 if w is on 72 ; and symmetrically 

Pj^{Sg{.,w)) = Sg{.,w), P^'{Sg{.,w) - Si{.,w)) = Sg{.,w) - Si{.,w) 

on 72 if w is on 71. The corresponding identities for Sg also hold. Together with the earlier convergence 
result for P^^ , s S {i,o}, we get for instance that 

\\pf^iRB{Sg{.,w)) - RB{Sg{.,wm^ = 0{fjiS)) 

(uniform norm on 7^), uniformly in w S 72. 
Consequently, if w is near 72, we have 

Sg{b,w) = ^Rb (e*^('")5g(6,w;)) + ^Rb [e-"''-^^ Sg{b,w)) + OifjiS)) 

for b a black vertex in a compact subset of A \ 71, and in particular for b on either side of 71 (ie both 
definitions of Sg agree up to 0(^(5)) near 71). It follows that KgSgl.jW) — is supported on 7^ and 

is 0{5fj{S)) there. Thus {KgSg), seen as an operator C"^ (where 7^ = U 7^) is such that 

\\\KgSg - Id|||ii(^«.) = 0{fi{S)) (since \-/^\ = 0{S-^)). Thus for S small enough, KgSg : O"' C^"' is 
invertible, and \\\{KgSg)~^ - Id = 0{ii{6)). Then for w € , we set 

S'g{.,W)= Yl {KgSg)-\w,w')Sg{.,w') 

SO that KgSg{.,w) = Syj. It follows that 

Sl{b,w) = ^Rb (e-(-)5,(6,«;)) + ^Rb (e--(-) 5,(6, ^i;)) + 0(^(5)) 

uniformly 'm w G 7^ and 6 in a compact subset of ^4 \ (71 U 72). 

For w in general position in A, w is in the exterior of 71 or in the interior of 72 (or both) . The two cases 
are similar, so assume that t« is in a compact set of A \ 71, on the; exterior of 71. Then we define Sg{.,'w) as 
we did for w near 72. Then KgSg{., w) — is supported on 7}^, is o{6) there, and we may set: 

S'g{.,w)=Sg{.,w)+ {K,{Sg{.,^m^')S%^o') 
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so that KgSg{.,w) = 6w Wc know that there is a most one / vanishing at infinity such that Kgf = 6^ 
(for S small enough). Thus Sg{.,w) does not depend on the choice of 71,72- By moving 71, 72 towards the 
boundary cycles of A, we can extend estimates of Sg{b,w) for {b,w) in a compact subset of \ Aa- 

Finally for w outside of A, it is easy to construct Sg{., w). For instance if w is in Do\Di, one starts with 
a truncation of K~^{., w) and set 

S'gi.,w) = l^...K-\.,w)- ^ Kgil^.^.K-\.,w))S'g{.,w') 

where 7'"* denotes the inside of the closed cycle 7. 

Let us summarise the results of this section. Recall that M is a bipartite graph derived from a rhombi 

tiling As with edge length S going to zero along some sequence; it is equipped with a linear operator 
K : C*^^ — >■ . We consider Di (resp. Dg) a simply connected neighbourhood of (resp. 00) in C, 

such that A = Di (1 Do is an annulus separating from 00. The (sequences of) graphs Sg = Sf , s G {i, o}, 
are bipartite graphs equipped with Ks : C^" — >■ C"" , a nearest neighbour (or finite range) linear operator. 
The pair {Ei,Ki) is obtained by modifying (M, K) in D^ \ Do, and vice versa for {Eo,Ko)- In particular 
(Sj, Ki) and (So, Ko) agree with (M, K) in A. The glued data {Eg, Kg) agrees with {Eg, Kg) in Dg, s G {i, o}. 
We also assume that for any w G Ef^ , there is a unique K~^{.,'w) G C~' vanishing at infinity such that 
Kg{K-\.,w)) = 6^, s G {i,o}. 

Lemma 8. Assume that Si,Si, So, Sg ■ A'^ \ A a C are such that 

K-\b,w) = ^RBie'''^'"^Sg{b,w)) + hlBie-^''^-'^Ssib,w)) + Oir,{6)) 

uniformly in compact sets of A^ \ Aa, s G {i,o}, with {z,w) 1— )■ Si{z,w) — ^^^^^^ (resp. Si{z. w) =L=^ 

in A? , and lim^x^o '7(^) = 0. Assume that there are Sg, Sg : A'^\ Aa — ^ C with same regularity such that: 
if "f is a fixed simple cycle in A disconnecting from 00, then for w G A outside 7, Sg{., w)\j is in the inner 
Cauchy data space Ci defined by Si and {Sg(.,w) — S'o(.,w))l-^ is in the outer Cauchy data space Co defined 
by So; and the corresponding conditions for Sg and for w G A inside 7 also hold. Assume that CiHCo = {0}, 

c^nC'o^io}. 

Then for S small enough, for each w G E^ there is a unique Kg^{.,w) G C"9 vanishing at infinity such 
that Kg{K~^{.,w)) = Syj, and 

K-\b,w) = ^i?B(e^^("')5,(6,«;)) + hlBie-'"^^^ Sg{b,w)) + 0{ri{5)y 5\\og6\) 

uniformly in compact sets of \ A a, s G {i, o}. 

Let us remark that in the case where Eo is actually bounded, the "vanishing at 00" condition is void. 
Conceivably, for unbounded Eo (agreeing with M far enough) , one might find use for "divisor" type boundary 
conditions at infinity: f{z) = 0{\z\"') as z — )• 00, n e Z fixed. In this case the previous surgery argument is 
still valid. 

7 Electric correlators 

We are now interested in scaling limits of vertex correlators. For instance in the plane, one may consider 
asymptotics of 
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where J2 = 0. Heuristically, we expect these asymptotics to be governed by electric correlators for a free 
field (: exp{iJ2j Oij4'{zj) ■)i with g = tt~^ (the value read eg from Corollary j4| . However, the discrete height 
function is at every point deterministic modulo Z (in the normalisation chosen here). Thus 

(a,) ^ {e'^^ 

is 27r-periodic in each variable, which is not the case of the scalar free field electric correlators: 

(: exp(2j7r^Sj-(/)(zj) :)c cx JJ |zj - Zjp"'''^" 

j i<3 

This may be seen as a manifestation of the compactified nature of the height field. 

The relevant Cauchy-Riemann operators are those associated to the line bundle Lp over the punctured 
sphere E = C \ {zi, . . . , where p : 7r(E) — )■ U is a unitary character. The two types of variations we shall 
consider are the isomonodromic family (zi, . . . , z„) Lp{zi, . . . , Zn), and Jacobian family p i— >■ Lp. 

The analysis relies on a rather precise description of the corresponding discrete operators and their 
inverting kernels, in particular near the diagonal. 

For notational simplicity we take 6 = 1. 
7.1 Discrete holomorphic functions with monodromy 

We proceed with a local study of discrete holomorphic functions and inverting kernels in the presence of a 
singularity. As before, we consider a rhombus tiling A of the plane. Let us mark the midpoint vq of an edge 
of A (that is, the center of a face of M). Up to scaling and centering we may assume S = 1, vq = 0. 




Figure 2: Local geometry of AI near the singularity vq (black disk) 

The unitary characters of 7ri(C\{tio}) — Z are identified with the unit circle U. Fix a non iriwaZ character 
X (identified to an element of U 7^ {1}). 

We consider (C*''^^)^ the functions on the lift of Mb to the universal cover of the punctured plane C\{'i;o} 
which belong to the character x- Elements of (C*^^)^ may also be seen as multiplicatively multivalued 
elements of C*^^ . The space (C*'^"')^ is defined similarly; we also consider the restrictions (C^^^)^ and 
^(^Mj^-j^^ We are concerned with the operators 

At ■■{C^^^)x ^ (C*-^^)x 
K :(C*^«)x (C*^»')x 
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We begin with a few basic estimates on harmonic functions. In what follows, B{0, R) = {z E My '■ \z\ < 

R}. 

Lemma 9. 1. There is e ~ £(x) > such that for n large enough, if f G (C*^^)^ is harmonic in -8(0, 2n), 



then 



sup |/(x)|<(l-e) sup 

reS(0,n) xedB{0,2n) 



Similarly, if f £ (C*^^)^ is harmonic in A{n, 3n) = B{0, 3n) \ B{0, n), then 



sup 

(|n,5, 



|/(a:)|<(l-e) sup \f{x)\ 

x^dA{n,3n) 



2. If f E (C*^^);j- is bounded and harmonic, / = 0. 

3. If y E My , there is at most one function G^(., y) G (C*^^),,- vanishing at infinity such that ZVGt^-(., y) — 
5y (here Sy designates an element of (C*^^);^ such that 6y{y) G {x" : n G Z} and 6y{x) = otherwise). 

4- If y & My, G^{.,y) exists and satisfies 



G^{x,y)^Oi 



) 



•if \x-y\> \y\/2, for some e = e{x) > 0. 
5. For all x,y E My, G^{x, y) = G^{y, x). 

6- ify-y', 



G^{x,y')-G^{x,y)^Oi\y\-\ 



)) 



if \x — y\ > |y|/4 and G^{x, y') — G'^{x, y) — 0{\x — y\ ^) otherwise. 



Proof. 1. Take x G dB{0,n); up to rotation, we may assume arg(x) = 0. Consider 

2 4 

C = -n) n {z : I arg(z)| < tt - Eq}, 

Eo > small enough and fixed. From the convergence of discrete harmonic measure (see [7], Theorem 
3.8), we may deduce that there exists 77 > such that for n large enough, the probability that the 
random walk on F started from x exits G on either the top or bottom side of the cone {z : \ arg(z)| < 
TT — £0} is at least rj; let us denote T and B these events. If ?/ is a point on the top side, a random walk 
starting from y disconnects the bottom side from dB{—x, |) before exiting B{—x, ^) with probability 
at least rj' > (uniformly in n large enough, y G j4(|n, this may also be seen using harmonic 
measure estimates). 

Hence with probability at least 77' > 0, we may couple the random walk started from x conditional on 
T with the random walk conditional on B in such a way that they couple before exiting i?(0, 2n) and 
their winding around differs by 27r. This shows that 

|/(x)| < (W|l + xl + (l-2w')) sup \f{x)\ 

xedB{0,2n) 

and since x 7^ Ij we have |1 + xl < 2. The same argument works in the annular case. 
2. By iterating 1., we get 

\f{x)\<{l-er sup |/(a:)|)<(l-£)"||/|U 

xedB{0,\x\2^) 



and consequently f{x) = 0. 
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3. Follows from 2. 

4. Up to rotation we may assume arg(y) = 0. Let Ce = dB{y, \y\/2) and Ci — dB{y, \y\/4:). For x € Ce, 
y G Ci, let /ie(t/, {a;}) be the harmonic measure on Ce seen from y € Ci and be the chiral 
harmonic measure on Ci seen from x € Ce- Explicitly, if C" (resp. ?/„) is the lift of Ci (resp. y) to the 
n-th sheet of the universal cover of C \ {0}, 

IJ-i{x,{y}) = X! x"'Harmuc^{x,{ym}) 
mez 

Reasoning as in 1., we can couple a random walk starting from x which approaches the crosscut (— oo, 0) 
from above before reaching C'i with a random walk approaching it from below in such a way that, with 
probability bounded away from 0, if the lift of the first random walk exits at y„i, the second one exits 
at j/m-i. This shows that ■)\\tv < 1 ^ £ for some constant e > 0. 

If C^{.,y) exists, it satisfies (say for x G B{y, |j/|/4)) 

Gxi^^y) = GB{y,\y\/2{x,y) + ^ He{x , {a}) n^{a, {b})G^{b, y) 

Denoting by T : L°°{Ci) -> L°°(C,) the operator: 

(T/)(x)= J2 Me(a;,{a})/i.(a,m)/(&) 

we have ||T|jioo_j.^oo < 1 — e since OIItv < 1 ~ e for all x. Conversely, we can set on C^: 

and extend it harmonically outside Ci; inside Ci we add the harmonic extension to Gb(j, |j^|/4(., y). It 
is easy to check that the resulting function satisfies the defining condition for G^. 

We have CBiy,R){x,y) = 0(1) for x e B{y,R)\Biy,R/2) and GB(j,,fl) (x, y) = 0(log|i?|) for \x-y\ « 
\y\ (see [7j). Thus G^{x,y) = 0(1) for x e d. Consequently, by 1. we have 



Gx(a;,y) = 0( 



) 



ii\x-y\ < \y\/2. 



5. By uniqueness 3., G^{x,y) is x-multivalued in y for x fixed (reasoning on the universal cover). To 
evaluate 

^C^ix,y) 

we fix y and consider it as a x-multivalued function in x. Then observe that it decays at infinity and 
has the same Laplacian as a lifted Dirac mass Sy{.). Consequently by 2., 

^Gx{x, y) = 5y{x) = (5^(2/) 

which concludes. 

6. We write G^{x, y') — G^{x, y) = G^{y', x) — G^{y, x). If |a; — y| > \y\/2, a combination of the estimate 
in 4. (for G^) and a discrete Harnack estimate concludes. If a; e B{y, |y|/4), as in 4. we may write 

GxiVi x) = Gr{y, ^) - ^ log l^^l + Hv) 

where Gr is the "free Green function" ([7 , Theorem 2.5) and h is harmonic in B{x, \x\/2) and uniformly 
bounded. We conclude with the Harnack estimate and asymptotics for first differences of Gr. 

□ 
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We now turn to discrete holomorphic functions. 
Lemma 10. 1. The space of x-iT^ultivalued bounded discrete holomorphic functions: 



{/e (C^^«)x : \\f\U<oo,Kf = 0} 
is one- dimensional and is spanned by a function f^ with the following asymptotic expansion: 

f^(u) - Rb (2^-ir(l - s){u - voy-^) + Rb (f2-^r(.s)(7i-t;o)"') + 0{\u - vo\''-^ + \u~ ^or^) 

where s g (0, 1), x = e^*^*, andr = (xo — vq)/\xo — vq\ if xq is the vertex ofT adjacent to the singularity 

Vq. 



2. More generally, for s G (0, 1), k (zN, there exists fk,x discrete holomorphic and x-i^^ultivalued s.t. 
fk,^{u) = Rb (t2^-'T{s - fc)(u-i;o)'~') + 0{\u - v^r""-' + \u- iK,\''-'-^) 
with fo_-^ — f^ and for b a black vertex adjacent to the singularity, 

[b — vo)^ 1 — e ^^'^^ sin(7rs) 

Remark that fk,x is also discrete holomorphic and ^-multivalued, and that = '''fx- ^s also point 
out that if is a black vertex adjacent to the singularity, 

ie plugging b in the leading term of the asymptotic expansion of fk,x yields the correct value up to the 
positive multiplicative constant i,\ ■ i — n- 

Proof. 1. Let / e (C*^^) be bounded with K/ — 0. Then a local computation shows that f\Mv is 
harmonic except at a;o, the vertex of F adjacent to the marked edge of A. Consequently, f\Mv decays 
at infinity and is proportional to G;^.(.,a;o). Besides f\MF is locally given as the harmonic conjugate 
of f\Mw and is thus given modulo a global additive constant (on the lift to the universal cover). This 
additive constant is determined by the fact that f\MF is x-^iultivalued. Thus the space of bounded 
harmonic functions in (C*^^) is at most one-dimensional. 

We need to exhibit a non trivial holomorphic function. Adapting the argument for the construction 
of in [55] (see also [3, Appendix A, which we follow more closely), we use an integral representa- 
tion involving the "special" discrete holomorphic functions (or discrete exponentials). Recall that the 
singularity is located at vq = 0, the midpoint of an edge of A which abuts xq (E T and yo G (thus 
Xq = —yo)- For a parameter A, we set 

exixa) = 



and for x E T , y E adjacent. 



1 - 2(^0 - ya) 

1 - 2 - xo) 

exjx) ^ 1-liy-x) 
exiy) l-^{x~y) 
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It may be checked that e\ € (C*'^^) is well-defined and discrete holomorphic in the sense that Ke\ = 
(recall that K is obtained from if by a gauge change). For u S Mb, we have e\{u) = 0(1/ |A|) as 
A — )■ oo. Besides, one can choose a path (uq . . . w„) on A fi'om uq = xq or yo to u„ = u in such a way 
that uq,ui — uq, . . . ,Un — Un_i lie in the cone: 

{z : I arg(z) - arg(M)| < tt - Eq} 

for some eo > 0. Consequently the poles of e\{u) are in {A : | arg(A) — arg(M)| < tt — Eq} (and have 
norm 6/2). 

For s e (0,1), set 

nUOD 

fs{u) = / A-^e_A(^^)dA 



"'0 

which is X = e^"'*-multivalued. Here /g"°° represents integration along a half-line from to infinity 
in the direction u. If w G Myi/, ui,...,U4 its four black neighbours, one can show that the same 
integration path may be used for ui, ... ,1*4 and consequently by linearity {Kfs){w) ~ 0. 

For small A we have 

exiu) = exp(AM + 0{\uX^\ + jA^])) 

and for large A 

exiu) = -— ^ . exp(4^ + Oi\u\-'\ + \X-'\) 

where e = 1 on F and e = — 1 on F^^ (as in [3S]). (Notice that yz| = exp{2e + 0{e'^))) . The asymptotic 

expansion comes from small values of A (say |A| < l/y^|it|) and large values (say |A| > ^/H). The 
intermediate values of A contribute to an exponentially small error in the asymptotic analysis. We 
have 

/•uoo 

/ A^" exp(-Au)<iA ^ F(l - sy-'^ 
Jo 

uoo / ^ 







A"" I : ^exp(-4uA-i) dX = r(s)(4u)- 

A(xo -vo) J xq-vq 



The error terms are estimated via 

t-'exp{^\u\t)0{\u\t^ +t')dt^ / t-'exp{-t)0{t'^ + ^0{\uY-^) 



JO m 

and similarly for the other integral. This shows in particular that fs is not identically zero. We then 
obtain by multiplying by a constant and applying a gauge transformation. 

2. For ken, se (0,1), set 

fk.sW) = I X''-'e-x{u)dX 

where 7 = 7(u) is an integration contour consisting of a large counterclockwise circle around zero 
(encircling all poles of the e.(M)'s); two connecting segments on t2(0, 00); and a small clockwise circle 
around (with all poles on its exterior). As before it may be checked that Kfk.s — 0, since for a given 
w € Mw, the same integration contour may be used for all black neighbours of w. 
We choose a lift of the mapping u i-> A = — to the covers oi {u : u ^ v^} and {A : A 7^ 0}. We 
choose determinations of log on these pointed covers in such a way that for A = u^^, log(A) = — log(u). 
Then 7(u) is a simple contour on the cover of {A : A 7^ 0}, for instance by specifying that the outward 
segment is in the direction u^^; then the inward segment is on the next sheet of the cover. 
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Notice that, up to rotating the lattice around vq, one may assume that it — f o G (0, oc), in which case 
one may a fixed branch of log, which will be convenient for local computations. More precisely, if 
M' = e'^M, one may choose branches in the definitions of in such a way that f^^ {e^^u) = 

In order to obtain an asymptotic expansion for large R = |u|, one may take the outer circle of the 
integration contour with radius R and the inner circle with radius R~^ . The contribution of the inner 
half of the contour 7 is, modulo exponentially small error terms, 

^,_fe_i I x^-''e-xp{-Xu/R + 0{uX^/R^))dX = 0{\u\'-^-^) 

where ji is a contour consisting of a clockwise unit circle around connected to infinity by two half-lines 
along u{0, 00). Similarly, the contribution of the outer half of 7 is: 

^{k-s) I ^k-s ( ^ exp(-4M(i?A)-i + 0{u{RX)-^ + {RX)-^)] dX 

J -To \X{xo-vo) J 

where 70 consists of the unit circle connected to zero by two rays along u{0, 00). 
Up to rotation and inversion, we have to evaluate 



/ 



where 70 is a contour consisting of a clockwise circle around connected to infinity by two half-lines 
along (0,00). (For definiteness, set log(re*^) = logr + i9 for S [0, 27r), and accordingly = 
exp((s' — 1) log(t))). We observe that, by contour deformation. 



/ 

J 'To 



t''-^e-*dt = (1 - e2^'^"')r(s') 



when s' G (0, 00); and that both sides are analytic in s' (the RHS has removable singularities on — N). 
Consequently they agree for all s' G C. Thus 

fkAn) = e ^^~^~'^ T{s - fc)(4w)'=-^ + 0(|«|'=-^-2 + lur'^-') 

Xq — Vq 

On the other hand, for u G {xq, yo}, the residue formula yields: 

fkA^)^ [ A"-^— ^7^^ .dA = W^('=-^)(u-z;o)^-'=-^ 

J l + X{u-vo) 



□ 

Lemma 11. Let w € Mw be one of the two white vertices adjacent to the singularity vq- Then there is 
Qw S (C-'^^);^ which is discrete holomorphic outside of w, has asymptotic expansion 

9w{u) = Rb (2'e"'^^^r{l - s){u - voY-') + Rb {2-'e-'''^^^T{l + s),!^"^-!) + 0(|u|-^-3 + \u\''-^) 

and such that 



Remark that for s = Q, g^j = 27rK (.,w). 
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Proof. We proceed as in the construction of /^^ , simply changing the normahsation of the discrete exponentials 
(or rather reverting to the normalisation used in [5S], [7]; we follow closely Appendix A of [7] here). Let 
xq e r, e r''' be the two black vertices adjacent to the singularity; let x'q € F, i/g e be such that w 
corresponds to the edges (xqx'q) and (yoJ/o)- ^"^^ A G C, set 

" ii~lixo-yomi-lixo-y'o)) 

and for a; e F, y e F^ adjacent, 

exjx) ^ l-^{y-x) 
exiy) l-^{x-y) 

Then Ke\ = and we have the following estimates: 

ex{u) = exp(A(u ~ w) + OdwA^I + jA^D) 



exiu) = . ^ ^ exp(4\^ + 0{\uX'^\ + \\-^\) 

{xo - yo){xa ~ y'^jX^ A 

For notational simplicity, we now assume that w — Q. 

Fix a ray ■i2o(0, oo) from w = to infinity which is disjoint of the rays m(0, oo) for all u € M, and such 
that arg(uo — wq) G (— tt, tt). We also fix a branch of log in C \ uop, oo) and set = exp(slog(z)) there. 
By rotating the lattice we may assume that mq = —1 and log is the usual branch in C \ (— cx),0] (so that 
{z~^y = z"" outside of the branch cut). 

Then set ^ 

~9{n) = I {-\)-'ex{u)d\ 

J —uoo 

which defines a single valued function on Mb- It can be identified with an element of (C*^^)^ by taking 
Mo(0,oo) as a branch cut. As before, one then checks that Kg vanishes except at w, by deformation of the 
integration path. 

We have 

i-Xy exp(Au)dA = F(l - s)u"^i 

DO 

i-X)-^ (- ^ — exp(4^A-i)) dX = ^ ^r(l + s){Au)-^~' 

\[xa-yo)(xo-yo)X^ ) [xq - yo)(xo - yo) 

Observe that (j/o — xo){yQ — xq) — e'^^"^^^ (if S = 1). We need to evaluate g at neighbors of w. In order 
to evaluate /"_^(— A)^''/(A)(iA, where / is rational of degree —2 without poles on the integration path, we 
may write (here x — e^"^) 



ix-' - 1) / (-A)-V(A)dA = i i-X)-sf{X)dX 



where 7 is a closed counterclockwise contour containing all poles of / and not crossing the ray — m(0, 00) 
and A 1-^ (— A)^-* is the determination of A M- (— A)^'' with a branch cut on — m(0, 00) which agrees with the 
reference determination on the right handside of —u{oo,0). Specifying to /(A) — (^i_xzi){i-\z2) ' residue 
formula yields 

J-uoo (1 - Azi)(l - AZ2) X - X ){Z2 - Zi) \ 
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Let us denote vq,vi,V2, V3 the midpoints of edges of the rhombus of A corresponding to w (where vq is the 
midpoint of (xoj/o)) hsted counterclockwise, and also (6oj ^2, ^3) = (a^Oi 2/0: a^oi 2/0)- We have 

1 

e\{bi 



(1-Aw,_i)(l-Aw,) 



for j = 0, . . . , 3, with cyclical indexing. For definiteness, let us assume that uq was chosen in the interior of 
the cone generated by xq and vq. Then we obtain (in terms of the reference determination) 

ii-vsY i-v2Y) 



9{bo) = 






2i7r 




(1- 


X' 






gih) = 






2i7r 




(1- 


X" 




~ «o) 


= 










(1- 


X 




-fi) 








2i7r 




(1- 


X 




- W2) 



Taking into account K(w, bj) — i{vj-i — Vj), we get 

{Kg){w) = 27Ti~voY 

In order to obtain g^, we multiply g by 2^e^^'^^^ and apply a gauge transformation. □ 

Lemma 12. Let x = e^''"' , s € (0, i). Let f £ (C*^^)^ &e bounded and s.t. K/ = owteide of B(0,R). 
Then there exists c — c{f) such that for eg > 0, [z] > i?, 



f{z)^cf^{z) + O{\\f\Bi0M)\\oo\z/R\ 



More precisely, for s e (0, \], there exists Ci(/) = 0{\\f\B(t)M)\\ooR'') , C2(/) = 0{\\f\B(o,R)\\ooR^ ") such that 
foreo > 0, \z\ > R, 

f{z) = c^f^iz) + 0(||/|B(o,fl)l|coN/i?r-^) 

= Ci/;,(z)+C2ffx(^)+O(||/|B(0,fl)||ook/i?r^-'+^") 

where = g^ for a w £ Mw adjacent to vq. 

Proof. Fix £, 77 > 0. Then there exists i?o > such that for all g discrete harmonic outside of B{0, Rq), there 
exists a continuous harmonic, x-multivalued function g s.t. for all z € My, \z\ > (1 + i])Ro 

\g~ g\{z) < e\\gioB{o,R.o)\\oo 



Moreover Rq is uniform in the choice of rhombus tiling A under (3.1). This is obtained by contradiction 
using equicontinuity and Ascoli-Arzela (see |7] for similar arguments). Similarly, if g is discrete holomorphic, 
then g\Mp converges to a continuous harmonic function g^ , which is conjugate to g. 



Let g be bounded, harmonic, ^-multivalued on {z : \z\ > 1}. Then g can be written locally as the sum of 
an holomorphic function h and an antiholomorphic function h. The decomposition g — h -\- h is unique up 
to additive constants. It is easy to see that the additive constant can be set so that h, h are x-multivalued. 
From the Laurent expansion of /i(z)z~*, h{z)z''~^, we get that 

g{z) = + E f>nz~'-'' 

n>l n>0 
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for coefficients ((x„), which are bounded by ||5||as(o,i)l|c!o- 

Thus for all 1^1 > C > 1, 

Let us now consider / (restricted to My for notational simplicity) . Let ^„ be continuous, harmonic such 
that for \z\ > (1 + r])Rn 

1/ - 9n\{z) < £||/3B(0,fl„)l|oo 

{e,Rn to be fixed later). There is c > 0, a„ = O(||sf|aB(0,R„)||oo-R*) such that for C > 1 large enough, 

\Z\ > CRn 

\gn{z) - anZ-'\ < \zr\l + cC-')\\giaB(0,il+r,)Rj\ < \Z/Rnr\l + CC-'){1 + e)\\f\OB^O,Rj\oo 

Given that fx{z) = Rb{z~^) + 0(|^|*~^) (up to normalization), we get that for z > CRn, 

\f - aMz) <{e+ \z/Rnr\l + cC-'){l + e){l + cRl'-^))\\fi9B(o,R„)\\oo 
Similarly, one can find coefficients /3„ such that 

1/ - Pn+lfM <{e+ \Z/Rnr\l + £'))!!(/ - /3n/x)|ai3(0,fl„)l|oo 

for l^l > Rn+i = CRn- Choose e small enough and C, i?o large enough so that e + C*~^(l + e') < c^~^~^'^° 
(and Ro> Rin any case). It follows that 

IK/ - /3n/x)|aB(0,fl„)||oo = Om\OB(OM\URn/Ror-'+'°) 

Since ||(/ - (infx)\dB{o,Rn+i)\\oo is of the same order, we get 

/3n+l -I3n = Om\aBiO,Ro)\\oo{Rn/R<,r-'+"'K) 

If we fix £0 > small enough such that 2s — 1 + eq < 0, this is summable. Moreover, if ^ = lim/3„, we have 

(3n = /3 + 0{\\{flOBiO,Ro)\\oo{Rn/Ror-'+"'R'n) 

and finally 

IK./ - /3,/x)|f?B(0,-R.„)l|oo = 0(|K,/|aB(0,iio)l|oo(^n/-Ro)*~^^'^'') 

In order to obtain the more precise statement, we write g{z) = boZ~^ + aiz^~^ + 0(|z|~*~^) we proceed 
similarly to show the existence of coefficients (/3n,7n) such that 

ll(/-/3n/x -7n5x)|aB(o,ii„)l|oo = 0(||(/|aB(o,fio) lloo(i?„/i?o)-^-'+^°) 

Thus 

(/3„+l - f3n)fx + iln+l - ln)gx = »( ||(/|aB(0,iio) II oo (i?n/i?o)-^- '+^° ) 

on 5(0, Rn+iY- This shows (by specialising at points close to the positive and negative half-lines respectively, 
say) that 

Pn+l -f3n = Om\OBiO,Ro)\\oo{Rn/Ro)-'-'^'°K) 

and subsequently 

7„+l - 7n = Oi\\if\aBiO,Ro)\\ooiRn/Ror'-'+'°Ri-n 

and one concludes as before. □ 
In what follows, we assume without loss of generality that s € (0, j). If s € (— g, 0), set = K^^. 
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Lemma 13 (Chiral Cauchy kernel). Let x = e^"*, s G (0, For w G Mw, there is a unique function 
K~^{.,w) e (C^^)x such that KK~^{.,w) = 5^ and 

for So > 0, \b\ > 2\w\. Moreover, K-^(6,u;) = 0{j^^{\w\/\b\y)) for \b\ < \w\/2. 

Proof. The uniqueness follows from the fact that there is no nonzero holomorphic function / s.t. f{z) = 
o(|z|-). 

If w corresponds to the oriented edge {xqXq) of F, we may set 

S{b) = G^{b,x'o)-G^{b,xo) 

for & € r. One may extend S to in such a way that S{yo) = {yo is the vertex of adjacent to the 
vertex singularity at vq = 0) and KS = except at w. We know that 

5(6) = 0(^((|6|/HrA(H/|6|r)) 

for 6e r, \b-w\> \w\/A. 

Consider a branch cut running from to infinity in the half-space {z : ^{zw) < 0}. Applying the discrete 
Green's formula in the slit plane (first intersecting it with a large disc) shows that the sum of (weighted) 
discrete derivatives of S across the slit vanishes. By harmonic conjugation, this is saying that S{b) as 
6 — >■ 00 on . 

By harmonic conjugation and the Harnack principle, we get the following estimates on S for |6 — w| > 1^1/4, 
b e r^, by integrating along a ray from yo to b if |&| < \w\ and from oo to |6| if |6| > \w\: 

sib) = oCj^mMrA{\w\/\b\r)) 

By the previous lemma, there exists /3 = 0{\\f\\QB{o,2\w\)\w\'') such that K^{.,w) = S{.) — /3f^ satisfies 

K^ib,w) = 0{\w\-\\b\/\w\y-') 
for \b\ > 2\w\, say, and K^(6,w) = 0{\w\-Ww\/\b\y) for |6| < |w;|, |6-w;| > \w\/4. 

□ 

Lemma 14. Let x = e^*''^ s e (0, i). Let f e (C^^)^ be bounded and s.t. K/ = in B{0,R). Then 

f{z) = O{\\fl0B(0,R)\\oo\z/RD 

Proof. Let / = /1s(o,k)- Then K/ = 0(||/|aB(o,ii)||oo) and the support of K/ is contained in white vertices 
adjacent to black vertices vertices belonging to dB{0, R) c Mb; let us use the same notation for this subset 
of Mw ■ Then consider 

/- E (K/)(HK^(«^,.) 

wedB{0,R) 

This function in (C*^^)^ is discrete holomorphic and 0(|2;|*~^) at infinity, hence vanishes identically. Con- 
sequently, we have the (chiral) Cauchy integral formula: 

m = E mM^xi^,b) = o{\\fieBio,R)\\oo\b/Rr) 

wedB{0,R) 

mB{0,R). □ 
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Lemma 15. Let x = e^^""' , s e (0, 1). 
1. The kernel satisfies: 



2 \'K{b-w) \w J J 2 \Tr{b-w) \w 
for < \b\ < 2\w\, \b-w\ > l\w\. 
2. If w is adjacent to the singularity vq = 0, we have 

and if m,oreover the edge [bw) is adjacent to the singularity, and vi is the symmetric image ofvo across 
(bw), then for b €T, Vi,b, bo in counterclockwise order around w, 



K{w,b)K-\b,w) 



1 - X ^ V V"' - ^0 



W — Vi 



Note that the conditions in 2. are not restrictive, by exchanging the roles of r,r''^ and/or conjugating 
the lattice. Also, remark that 

I _ gis(arg(u)-iii)-arg(ju-i;o)) ^ 

lim =2^;;^ = ;^(arg(u; - vq) - arg(w - vi)) = K(w,6)K" {b,w) 

Before proceeding with the proof, let us sketch an alternative construction for K^^{-,w), for w adjacent to 
Vq, from which one can deduce the evaluation K~^(6, w) for b, w adjacent to vq. Consider a branch cut 7 (a 
simple path on Af^) running from infinity to and vi. 7 = (. . . , v^2^ ^-i, ^^Oj ''^i)- Wc have constructed basic 
discrete holomorphic functions f-^^y^ with monodromy at Vi, i = 0, 1. Using the branch cut 7, we can look at 
them as single valued on the same fundamental domain. Then „j , seen as a function with singularity at 
Wo, is discrete holomorphic except at w. From the asymptotic expansions, we can find a (non trivial) linear 
combination of fx,VQ-, fx,vi which is o(|2;|~'*) as z — >■ 00. Thus, by the characterisation of K^^{.,w), it can 
be written as a linear combination of fx,vo^ fx,vi- Prom the exact evaluation of these functions near their 
singularity, one recovers 2. 

Proof. 1. Let us fix w e Mw and set R= \w\. We first construct an approximate kernel S G (C'^''^^)^. If 
b e S(0,ro), set S{b) = cfx{b). If 6 e B{w,ri), set 

S{b) = K-\b, w) + ^Rb(- -)+^Rb{ ^ 

2 V"" w J 2 V"" ^/ 

and otherwise 

^' 2 ''\TT{b-w)\wJ J 2 ''\7:{b-w)\w 
Here c, ro,ri are parameters yet to be specified. Then KS{.,w) = dyj in B{0,ro) U B{w,ri). Set c s.t. 

Then the discontinuity of S on dB{0, rg) is of order 0{rQ/\w\^~^^ + rQ~^ /\w\^^''). 

The discontinuity of S on dB{w,ri) is of order 0(ri/|wp). Elsewhere KS{.,w) is controlled by the 
third derivative of the continuous limit. 
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Next we use the a priori estimate on the kernel K^"'^ to correct S. First let us observe that K^^Cw, .) = 
S-K~^{KS). For this need only check that (KJ^^(K5))(6) = 0(|6|""^) as 6 ^ oo. The contribution 
from {KS) in B{0,2R) is easily seen to be 0(|6|'*~^). The other terms are handled as below . 
Let us estimate the correction K~^{KS) for b G A{R/2,2R), b ^ B{w,R/2) (where R = \w\, say). We 
simply add upper bounds (up to multiplicative constants) stemming from (KS){w') for respectively: 
w' € B{b,R/4); w' G dB{w,ri)- w' - w € A(ri,i?/4); w' £ A{R/2,2R) \ {B{b, R/A) U B{w, R/4)); 
w' i B{0,2R); w' £ 9B(0,ro); and w' £ A{ro,R/2): 

T.kR^ + j^-n + T.v^-n + Wji'-^'-li^-}^^''^R^ + ^''''^ )^ 

k=l k=ri k=2B. " 

Thus, up to a constant, we get the upper bound 



rf 1 To 
R^^r[R^R^ 



with ri = 0(Vi?), ro = 0(1). 



2. Here ^(., w) is proportional to gw Since (5 = 1, we have (dq — w)*(vo — w) =2 



□ 



We are now in position to state the analogue of Lemma 12 (concerning the asymptotic expansion of a 
discrete x-multivalued holomorphic function in a neighbourhood of oo) for the singularity at zero. 

Lemma 16. Let x = e'^"", s £ (0, 5). Let f £ (C*^^)^ be bounded and s.t. K/ = m B{0,R) so that in 
Ai^,R), _ 

f{z) = Rb{Mz)) + Rb{Mz)) + e{z) 
where i/'i(z) = Z]„>o i'^iz) = Z],i>o ^^s"- > 0. for z = 0(1), 

^ ?2^I>)-^^(") + r2-T(-.) ^^^"^ + 0(i?^+''||e|L,^(H,«) + i?^^-^+''||V.IL,A(f 



w/iere h^{z) = fi,x{z), and WiJjiW = HV'iH + ||V'2||- 
Proof. Let (p £ (C*^^)^ be such that K0 = in 5(0, r), 

(/)(z) ^ af^{z) + f3hx{z) + Rb(Mz)) +Rb{Mz)) + e{z) 

in A(|,r) = B{0,r) \ i3(0, |), say. Here (pi (resp. (/)2) is holomorphic (resp. antiholomorphic) and 
multivalued, with Laurent expansions (/)i(z) = z^ J2n>i "^"-^"i '^'2(2) = ^ * X]n>i ^n^"- Observe that 

^ ^ -J ^ I / z \ « (/>i(w)dw 1 / / z Y (l)i{w)zdw 

for z £ D{0, |r), and thus ||</'i||oo,c(o,eor-) < c(£o)||0i||oo,c(o,r), where c(£o) ^ ^o^^ goes to zero as Eq goes to 
zero (eq > small, to be fixed later). We have a similar estimate for ^2: ||02||oo c(o eoi-) — c'(eo)||</'2||oo c(o r); 
c'(£o)-er'- 

If 7r is a simple cycle on F approximating C(0, |r), let _B be the set of black vertices on or inside 75. We 
have the replication identity 

^{z)^ K{cP1b){w)K-\z,w) 
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where 7^ are white vertices which are on 7^ or are outside of 7^ and adjacent to a black vertex on 7^. For 

z € i?(0,r/2), say, we liave: 

4>{,z) = af^{z) + ph^iz) + 0(||£||oo,A(5,r)(|^|/r)-^) + ^ K((i?B(/.i + Rb<I>2)1b){w)K-\z,w) 
As in the surgery section, for z G ^(^,£0'^)) we may estimate the discrete Cauchy integral by: 



Consequently in A{^,sor), we have 

(/.(z) = af^iz) + ph^{z) + Rb{Mz)) +Rb{4>2{z)) + e{z) 

with 

for i = 1, 2, and 



(l|e||oo,A(5,r) + \\(l>i\\oo,A{i,r)r'^'' ^) 

Fix ?7 > arbitrarily small. For small enough eq, k > 1, we get 



and 



^+''11^11 < IIpII -4-IU.II ,r2s-l/-feN2s-l+l-s-») 



-0 "-"oo,A(iap:,eS+V) - "'"oo,AC-^,4r) 

SO that for r} small enough we get 

Ik 



^^eg. , , < {4)-'-" (lkl|oo,A(S,r) +C7.||<^,|U,^(;,.)r2-l) 

and for z = 0{1): 

cj>{z) = af^{z)+ph^{z) + 0{r'+^\\e\\^,Air^,r)+r^'^''-^\m^^ 

We have 

f^{z) = RB{f2-^ns)-z-^) + 0{\zr') 
h^{z) = Rb{t2'T{-s)z') + 0{\z\-'-^) 

Taking / as in the statement of the lemma, set (l)i{z) = 'ipi{z) — aoz^, (l>2{z) = tp2{z) — boz~^, a = f2-"r{s) ' 
^=V2^)- We get 

f{z) = af^iz) + ph^{z) + Rb{Mz)) + Rb{Mz)) + s'(^) 
where \\(l>i\\oo,A{§,R)^ and pR^ are of order \\il)i\\^^A(%,R)) and 

lk'lloo,A(f = Om\ooAii,R) + IIV'illoo,A(f 

which concludes. □ 

This applies in particular to / = K~'^(.,ii;), R = \w\/2. For b = 0(1), this leads to an estimate of 
K^^{b,w) (which is a priori of order 0{R'^^^)) within 0(i?'^''^^+''). By replication, for s < j, if (/) is discrete 
holomorphic and x-multivalued in a ball -6(0, R), we can estimate (/)(6) for b adjacent to the singularity based 
on values taken by (/> in the neighbourhood of a closed cycle around in A{^,R). 
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7.2 Pair correlations 

Let us now consider the situation where two points x (midpoints of edges of A) are marked. We consider 
functions with monodromy x ~ e^^^^, s € (0, |), around x and x around y. Wc arc interested in estimating 
Sp, a kernel inverting K operating on these functions, in particular near the diagonal. 

More specifically, we consider a family of rhombus tilings (A^)^ with respective edge length 6, 6 \ along 
some sequence. Two points xs,ys (midpoints of edges of A) are marked; we assume xg x, ys ^ y. Let 
p : 7ri(C\ {.T, y}) — > U be the unitary character s.t. p{^x) = X: p{lv) = X where 7^ (resp. 7^,) is a simple loop 
winding around x (resp. y) counterclockwise. Let {C^'^)p (resp. (C^^)p) be the space of multiplicatively 
multi- valued functions on Mb transforming according to p. 

First we consider the discrete Laplacian operating on My. 
Lemma 17. We assume x ^ y, p non trivial. 

1. There is no nonzero hounded harmonic function in (C^'^)p. 

2. For w € My, there is a unique hounded function Gp{.,w) S (C^^)p s.t. ArGp{.,w) — 5w. 

3. There a unique hounded continuous p-multivalued function Gpi^.^w^ on C \ s.t. ^Gp(^..,w^ — (5^ 
and Gp{z, w) ^ as z ^ x,y. 

4- As 5 \ 0, Gp{.,.) converges uniformly to Gp{.,.) uniformly on compact suhsets of {{z,w) e (C \ 
{x,y}f,zi=w}. 

Proof. 1. Let /i be a bounded harmonic function in (C^^)p. Fix z e My- Consider the random walk 
(X„) started from z stopped at r, the time of first return to z or the first exit of B(0,n), whichever 
comes first. Let Wx (resp. Wy) be the winding of {Xq, . . . ,Xt) around x (resp. y). An excursion 
decomposition of the random walk between successive visits to z before time r shows that: 

'^(^^ = 1-F(x/="X---H^.=^) ^^^^^-^'^- ^ 

As n — >■ 00, |E(x™==~"'!' = z)\ stays bounded away from 1 while V'{Xt = z) ^ 1. Consequently 
h{z) = 0. 

2. Uniqueness follows from 1. For existence, we set h{z) = E^(p(7)), where 7 is the concatenation of a 

random walk started from z stopped when it first hits w and a fixed path from w to z (one may also 
reason on a branched cover of My, stopping the walk when it hits a lift of w). Then h is bounded 
and harmonic except at w. It is indeed not harmonic at w by 1. Consequently it is proportional to 
Gp{.,w). 

3. By compactifying at infinity, a bounded harmonic function has a removable singularity at infinity. By 
mapping {x, y} to {0, 00} by a homography and lifting to the universal cover of C \ {0} via log, the 
problem it to find h harmonic except on w' + 2mZ s.t. h{w' + 2iiTk) = xh{w'). It is then easy to see 
that the solution to this problem exists and is unique (up to normalisation) and decays as ^{z) — )• ±00 
(corresponding to x, y in the original problem). 

4. Fix e > 0; we consider {{z,w) € {B{0,e~^) \ {B{x,3e) U B{y,3s)))^ : \z - w\ > e}. We need to show 
uniform convergence of Gp on this set. 

Let 7e, 7i be simple paths on Ts at distance 0{S) of the circles C{w,2e), G{w,e) respectively. Let 
B be the connected component of w in \ 7e. Let Tg (resp. r^) be the time of first exit of F^ \ 7e 
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(rcsp. \ 7i) by a random walk on F. Wc define (twisted) Poisson operators Pg : C'*'* C^* and 
p. : C7. ^ C^e by {PJ){z) = E,(/(X,J) and 

(P,/)(^)=E,(x^/(X,J), 

where A'' is the algebraic number of crossings of a branch cut between x and y (not intersecting 7e) 
by the random walk before Tj. (One may also reason on a branched cover of My, lifting / to a p- 
multivalued function). We omit the dependence on 6 for lightness of notation. 

As before, one can show that {PePi) is a strict contraction on L°°{ji), uniformly in w, 5 small enough 
(e > is fixed). Moreover, by starting a random walk on 7j and stopping it on its first return to 7, 
after its first visit to 7e, we get the following identity in L°°{'-fi): 

G%,w)i^, = G%{.,w)i,, + {P,Pi)G%w)\^, 

Consequently, 

00 

G%,w)\',. = C£^PePit)G'B{;W)\,., 
k=0 

a summable series in L^{ji). 

Clearly, the analogue decomposition holds for the continuous Green kernel Gp. We know that Gg{., w)^^. 
converges uniformly to Gb{-,w)\~^^ as (5 \ and in particular is uniformly bounded. Since (PePi) is 
uniformly strictly contracting, we only need to establish uniform convergence of each term in the 
expansion. 

We note that P^h converges to P,,h in L°° {^fi{uj)) uniformly in h. w as (5 \, 0, where is (the restriction 
to 7e(ui)) of) a 1-Lipschitz function. This may be shown by contradiction, using equicontinuity (of h\^. 
and of Pg /i on compact subsets of B). 

Similarly, Pfh converges to Pj/i in L°°(7e(to)) uniformly in h,w, where h is (the restriction to 7e(w)) 
of) a 1-Lipschitz function. We reason as above, the contradiction coming from the uniqueness in the 
following Dirichlet problem: find ho harmonic, bounded, p-multivalued in C \ {{x,y} U C{w,e)), with 
boundary condition a continuous function h on C(0,£). As in 3., this may be recast as a standard 
Dirichlet problem in C. 

Using iteratively these convergence statements for P^ , P/, we obtain termwise convergence in the series 
expansion for above. 

□ 

As was the case with one singularity, we can use to construct an inverting kernel Sp for K operating 
on {C^^)p. Among such kernels, it is characterised by the fact that Sp{.,w) restricted to My is harmonic 
except at endpoints of the edge of T corresponding to w. The restriction of Sp to My is Gp(., x') — Gp(., x), 
where w G Mw corresponds to the oriented edge {xx') of F^. On Mp, it is deduced by harmonic conjugation; 
it vanishes on the vertices of F^ corresponding to x, y respectively. A discrete Green's formula argument 
(cutting the domain along a branch cut from x to y) shows that these two conditions for Sp{.,w)\MF ars 
consistent. Indeed, the variation of the harmonic conjugate between these two vertices is proportional to 
the flux of Sp across a branch cut connecting the vertices; as Sp decays at infinity, this flux is the integral 
of the Laplacian of Sp in the fundamental domain defined by the cut, hence zero. Note also that Sp is not 
invariant under duality F -H- F^ (ie reasoning on the random walk on F^ produces iSj, ^ Sp). 

In order to estimate Sp, let us evaluate d^Gp, d^Gp (in the continuous limit). As a function of z, d^Gp is 
locally the sum of a meromorphic and an antimeromorphic component; the decomposition is unique up to an 
additive constant. The additive constant can be specified uniquely in order to make these two components 
individually p-multivalued. Moreover, Gp{z, w) — log \z — w\ is harmonic in 2; in a neighbourhood of w. If 
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^ _ g2j7rs^ g g -|^^^ ^Yiis leaves the only possibility: 

1-s 



dwGp{z,w) 



z — X \ I z — y 



w — xj \w — y J Tr{z — w) 



as there is no (nonvanishing) /9-multivalued, bounded holomorphic (resp. antiholomorphic) function on 
C\{x,y}. 



Wc note that Gp{z, w) = G^{w, z) and thus we also obtain convergence of discrete derivatives of Gp{z, w) 
in cither variable. 

We are seeking an inverting kernel Sp which vanishes at infinity: Sp{b,w) — )■ as 6 — )■ oo. Such a kernel 
differs from Sp by a discrete holomorphic function (freezing the second variable). 

Observing that the restriction to My of a /^-multivalued discrete holomorphic function is harmonic except 
possibly at the two vertices adjacent to x,y, we conclude that the space of bounded p-multivalued discrete 
holomorphic functions is at most two-dimensional. Let us construct such functions. 

Let 75 be a simple cycle on F at distance 0{6) of 7 = C{x, \x — y\/2). Let 7] be the outer boundary of 
the union of faces of corresponding to vertices of F on 75. Set f{z) = 2^r{s)~^ f^^x{z) for z G Mb inside 
or on 75 and f{z) = otherwise. Then: 

gAz) = m- Yl mMSpiz,w) 

weMw 

is p-multivalued and discrete holomorphic. Moreover K/ is supported on vertices of Mw corresponding 
to edges of 75 and 7^. Assume here s G (|, 1) (so that the dominant terms are holomorphic rather than 
antiholomorphic; otherwise conjugate). For z G My, we have 

Y,{mMSp{z,w)= I i{w-xy-^d^Gp{z,w) + o{l) 
{\^f){w)Sp{z, w) = J^{w- xy-^ * d^Gpiz, w) + 0(1) 

Taking into account {w — xy~^ = (w — xy~''r'^'^~'^ on C{x,r), we may write these integrals as contour 
integrals of closed forms and deform them to 7 to obtain 



- ^ (- ^r-' ( ( ^ V f ^ V"^ zT^^I + 0(1) 



g:o{z) = {Z- xy '^X\z~x\<\x~y\/2 

' z — X^ 

— X J \w — y) 'k{z — w) 

+{z-xy-^X\z-x\<\x-v\/2)+0{l) 

by the residue formula. Symmetrically, one may construct another discrete holomorphic function with: 

-1-s 



gy{z) = {Y^y''(^^^ +0(1) 



The asymptotics on Mp are obtained by harmonic conjugation (or by a similar argument), the additive 
constant being fixed by the condition on p-multivaluedness. Convergence is uniform on compact subsets 
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away from x,y,j. By writing the Cauchy formula on another circle 7', one obtains uniform convergence on 
compact subsets of C \ {x, y}. 

Thus one may add a linear combination of g^, Qy to Sp{.,'w) to obtain a kernel vanishing at infinity. 
Indeed, a function on My (resp. Mp) which is bounded and (discrete) harmonic in a neighbourhood of 
infinity has a limit at infinity, as is easily shown eg by a random walk coupling argument. Convergence on 
compact sets is enough to ensure that 

lim g^{z) = {x~yy-^ +o{l) 

Si milarly , \vcsiz^oo,zeMF 9x{z) = i{x-yY~^+o(l), Yivci^^oo^zeMv 9y{z) = {v - xY^^+o{l), lim^_>oo,2eMj. 9v{z) 
—i{y — xY^^ + 0(1). Thus, for S small enough, one may find a{w), b(w) s.t. 

Sp{z,w) = Sp{z,w) - a{w)gx{z) ~ b{w)gy{z) 

vanishes at infinity. Combining with previous estimates on Sp{., w), gx, gy, we identify a(w), b{w) up to o(l). 
Let us summarise the previous discussion. 

Lemma 18. For s E (— S small enough, for w M^y there is a unique p-multivalued function Sp{.,w) 
vanishing at infinity s.t. KSp{.,w) = 6^. Moreover, 

Sp{z,w)^-RBii -) \]+7,Rb[[ -] +o{\x^y\-^) 

2 \\w-x z-yj ■k{z-w)J 2 \\w - x z-yj tt{z - w) j 

with uniform convergence on compact subsets of {{z,w) G (C \ {x,y}Y ■ z 7^ w}. 

The continuous holomorphic and antiholomorphic kernels are specified by the fact that they vanish as 
z — >■ 00, and have singularities of order at most |s| at x,y. 

Let us consider a more general situation. Let xi, . . . , a;„, j/i, . . . , y„ be marked points, si, . . . , s„ G (0, 1). 
We consider p : 7ri(C \ {xi, . . . ,?/„}) — t- U the unitary character such that pi'jxj) = Pilyj)^^ — Xj = e^'^'^"^ 
where 7^ is a counterclockwise loop around z G {xi, . . . , ?/„} with no other marked points in its interior. We 
are interested in the associated operator (C*^^)p — > (C^^"')p. 

Let us address the case n = 2 using the surgery formalism (Lemma [s]). Let 7 be a simple loop with 
{x\,yi} in its interior Ui and {a;2,2/2} in its exterior Uq. Let pi(resp. p2) be the character corresponding to 
weights (si,0) (resp. (0,52))- We have 



2 ^\\w-Xj z-yjj tt{z-w)J 2 ^lV"'-Xj z - yj J tt{z - w) 



where all pairwise distances are of order f . The Cauchy data spaces corresponding to the limiting continuous 
kernels are easy to identify: 

C^ = {/|, : f{z) = g{z) ■ ^^i^y , (5,),^^ = 0} n Up{^) 

\w-xi z-yij 

= {/|, : f{z) = g{z) (^—^ ■ , {gz)^u. = 0} H Lip{^) 

\w~xi z-yij 

Co - {/|, : f[z) = g[z) ■ ^^i^y , (g-^) ^ 0,5(00) = 0} n Ltp{^) 

\W-X2 z-y2j 

Co = {/i, : fiz) = g{z) ■ , (<7.)|c,„ = 0,5(00) = 0} n Up{j) 
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Note that in general Cj 7^ Cj. We check that C, fl Co = Cj fl Co = {0} and find the glued Cauchy kernels: 



n{z-w) j-j-^\w-xj z-VjJ 



Z — Xj w — 



niz-w) fj-^ \w-Xj z- Vj 
By induction on n (using surgery at each step to add a pair of insertions), we obtain 



2 y7^{z-w) j-J-^\w-Xj z-yjj J 2 yiriz - w) fj[\w - xj z - yj I I 

with uniform convergence on compact subsets of \ As, S = C \ {xi, . . . , Through surgery, we also 
retain invertibility (ie Sp is uniquely characterised by KSp{.,w) = S^, Sp{.,w) vanishes at infinity), at least 
for small enough S. 

We would like to gain an explicit error estimate for Sp. The following is a simple (and rather crude) 
estimate, which will be enough for our purposes. 

Lemma 19. We have 



' ^ 2 yn{z - w) f}^\w - X, z-yjJ j 2 y^iz - w) f},\w - Xj z-yj j ^ > 
uniformly on compact sets of \ As, S = C \ {xi,. . . , 

Proof. First we estimate Sp{h,w) for w near a singularity, 5 in a compact set of E. From 

K-\h,w)=0{\wr) 
for X = e^"^*, s G (0, |) (with singularity at 0), and the replication identity: 

we get the estimate: Sp{b, w) = 0{\w\^''^^^^) for w close to xi, 5 in a compact subset of S. 

Now fix Wo away from the singularities; we wish to estimate Sp{., wq). Let us set Sp{., w) = K^^(., wq) + 
|i?B(a) + 5-Rb(/3) in B{wo,r); and 

2 yr{z-wo)^J-^\wo-Xj z-yj J J 2 yTT{z - wq) fJ{\wo - Xj z - yj 

elsewhere (a,/3,r to be fixed). We have 

Sp{z, w) - Sp{z, w) = ^ K(5p(., wo))Sp{., Wo) 

On dB{0, r), we have K((5p(., wq)) = 0{5r'^) for an appropriate choice of a, (3. For r < |w— wqI << !> we have 

K{{Sp{., wo)){w) = 0((5^|w— wo|~^)- For w in a compact subset of S\{t«o}, we have K((S'p(., wo)){w) = 0(6"^). 
For \w — wo\ ^ I, K((S'p(., wo))(w^) = 0{S'^\w — wo|~^) (ensuring convergence of the RHS in the replication 
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identity). Finally, for w close to Xj, K{{Sp{.,'Wo)){w) — 0{S'^\'w — Xj\ I'^jI ^) (and similarly near yj). This 
leads to: 

oo 

Sp{z, w) ~ Sp{z, w) = 0{{r6-^){6r^)r-^) + 0(r '(54) + 0( X! k{^^{kS)-^){kS)-^) 

n <5-i 

+ ^ fc(54(A:,5)-l^^l-3(fcj)-l«.l) 

Thus setting r — -\/(5, we get 

5p(z,u;)~5p(z,w) = 0(<5i-2--'^^l^^l) 
as stated. □ 

We will need to analyse Sp{., w) when w is adjacent to a singularity. This is obtained by a combination 
of Lemma 16 and the estimates we obtained for Sp in Y? \ A^- Let us consider the "Robin kernels": 



z^w \ TTiz — W) -^-^ \W — Xi Z ~ Vi I TTlZ — W) I TT ^ — ' 



w 1 7r(z — w) ^^^\w — Xj z ~ Uj J 7r(z — w) j tt -^-^ \ w — Xj w — yj 



7r(z-w) -^^^^ Vw^-a;j 7r{z - w) I tt^^\w-Xj w - y^ 



(so that rp{w) = rp{'w)). 

Lemma 20. Ifb,w are adjacent to the singularity x £ {xi, . . . , j/„}, s = ±sj. the corresponding exponent 
(s = sj. for X = Xk and s = — s^, for x = yk), and 6 € F, we have 

Sp{b,w)-K^,{h,w) = ^ (^) (*e"'("')z) + 0(|x - 2/ri'^) 

where e — e(s) is positive for s small enough and 

z = ^ +y f = hm ( nrpiw) - ^ 

Proof. Let w be adjacent to the singularity x = Xfe (which is in a face of M); by translating we may assume 
Xfe = 0. Let r be of order 1 and small enough so that other singularities are outside of i?(0, ?■). Also set 
5 = Sfc, X = e'"^. We consider /(6) = Sp{b,w) — K~^(6, z/;)lB(o,r) which is discrete holomorphic except near 
dB{0, r) and is p-multivalued. The goal is to estimate f{b) for b = 0(6) (in particular for 6 adjacent to the 
singularity). We have the Cauchy integral formula 

f{b) = J2Spiw',b)K{f){w')= J2 Sp{w',b)K{K-\.,w)lBio,r))iw') 
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for b E B{0,r). From Lemma 15 and estimates on Sp, we obtain for b € A{j, |): 

9 TTlr — 7/n\-LJ. 



Tt{z — w) y ^-fj \W — Xj Z — Uj J \W — Xk 



r(l + .s)e-"-W— 1 / -Q ^ ^ - - ^ ' ^ z - a;^ 



2 7r(z - \ "Uj ~ z-yjj \w - Xk 

where = 2:'E„>o«»^"' V'2(^) = z"'' X;«>o ^"^"^ 

«„ . E(i_i)|!!:! f-^ + E (^^ ^ ) I (1 + om 

Sttw" I Xk-yk Xk-yjJ ' 



r(i + s)e-*'^("') 



2™ ^ \ Xfe - ?/fc ^ \xk - Xj Xk - yj 



and by Lemma [TBI we conclude: 



f{b) = —^^ f fb) + -h^{z) + 0(5l-4max, |s,|-,n 



for b = 0(1), 77 > fixed. Taking into account (Lemma [10|), we have for b adjacent to the singularity: 

x-i 



/x(6) = -—^2(^-6) 



hxib) = h.xib) = - 7-1 (^(^ - 
and r — |(5 — x) if & e F adjacent to a;. □ 



7.3 Variational analysis 

The previous local estimates may be used to estimate 

for x,y E at large distance, where h is the height function associated with a matching m of M and 
^ — g2iTrs jg gxed (s G (0, |)). More precisely, we are interested in estimating 
within o{\y' - y\/\x - y\). 

Fix X G and a simple path 7' from x to y' on M^^; the penultimate vertex of 7' is y: j' = {x, . . . , y, y'). 
We also denote 7 = (a;, . . . , y), the path stopped at y. Let Er = Er{-f) be the set of edges of A/ crossed by 
7, with the black vertex on the right hand side of 7; and Ei = Ei{'^) the edges crossed by 7, with the black 
vertex on the left hand side. Then (f33') 

h{y) - h{x) = ^uiind{j) + ^ leem - E "'■'^em 
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In order to define the winding number of 7, we draw it as a succession of smooth arcs connecting midpoints 
of segments (bb') (ie centers of faces of A/), where b E T and b' E , in such a way that 7 crosses these 
segments normally. 

Let : Mb — >■ be the operator defined by K^{b,w) — x^{b,w) if {bw) G Er{j); Ky{b,w) = xK(6, w) 
if (few) e Eiij); and Kj{b,w) = K{b,w) otherwise. 

Let us remark that if 7 is another path from x to y, then K;^ is conjugate to by diagonal matrices involving 
the winding number of the loop obtained by concatenating 7 with 7 (in reverse orientation) . Moreover, by 
using 7 as a branch cut (lifting functions on M to functions in (C*'^)p, where p is the unitary character of 
7ri(C\ {x,y}) defined as before), one can identify with K operating on (C*^^)p. Correspondingly, we get 
an inverting kernel : C*^"' C^^^ . 

Lemma 21. Kyi<r^ : Mw Mw is a finite rank perturbation of the identity and 

Proof. See the proof of Lemma [2j □ 

We now analyse the effect of taking the last step along 7'. Let {bw) be the edge of M separating 
y, y' G . Recall that Sj{., w) can be characterized as the only function in C*^^ vanishing at infinity such 
that K^S^{.,w) = Sw] S^'{.,w) is characterised similarly. Since K^i — is supported on w and 

{K,Sy{.,wmw) = ! + (!- x^'mw,b)Sy{b,w) 

(here — X (resp. x~^) if b is to the left (resp. to the right) of 7'), we deduce 

Sy{.,w) = (1 + {l-X^'mw,b)Sy{b,w))S,{.,w) 

In order to evaluate (^Hy')-H-^)) / (^^Hy)-h{^))^ we have to evaluate det(i^^/ K"^)/ det(is:^K"^). Let us 
observe that K-^/S-y — Id+(if^' — K^)S^ is a rank 1 perturbation of the identity and that 



as bounded operators in L?'{My/, ^) where p.{w} = (1 + for some £ > 0. This follows from SjK^ = Id 

(indeed, {S^K^)Sh — Sb is in the kernel of and vanishes at infinity, hence is identically zero). Consequently, 

det(KyK^ ) ^ jg^^j^^,^.^) ^ det(Id+(i^y-A'^)S'^) = l+{KY~K^){w,b)S^{b,w) = l+{x^^-l)K{w,b)S^{b,w) 
det(i4r^K ) 



since K^' — has rank 1. Then we use the key estimate (Lemma 20) 



■-7 -"-7 

£.2 



5,(6, w) - Kri^(6, - ^ (^) (*e-(-)/(y - x)) + 0{\x - y|-i-) 



together with ±zK(w, b) = y' — y, ^^^-^ = and the exact result 



K(u;,6)Kr],(6,^) = ^— I 1 



— f 1- 

-x\ V y - w 



— s > 



to obtain 

det(ii'yK"^) fy'-w\^'f^ ^2^fy' 



det(i^-yK ^) \y-wj \ \y ^ x 

Thus we can state 
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Proposition 22. For s small enough, there exists a c — c(A, s) > such that 



Proof. From the previous computation we get 

(^^h(y')-h(x)^ . \y-x\' 



c\x - y\ 



l + 0{\x-y\-'-n 



l^.^h(y)-h(x)\^ . \yl _ 2;|-2s2 

with £ > for s smah enough. Thus if we fix x and let y go to infinity, we get 

\l^/^(y)-H-))\^c{x,K,s)\x~y\-^^" 
Since \l^-xt^iy)-H^')) / l^-xi<v)-H^)) \ = 1 + 0{\x - x'\l\x - y\), we conclude that c does not depend on x. □ 

More generally, if {zi,...,Z2n} = {xiT--,yn\ are marked points, Si,...,S2„ > are small enough 
exponents with si + S2 = • • • = S2n-i + S2n — 0, we have (by "integrating" Lemma 201 



2n 



(exp(2i7r^ Sj/i(zj))) '-cJll 



\2siSj 



3 = 1 



in agreement with the electric vertex correlator heuristic. The pairing of insertions (zi, 22),- ■ • , (^2n-i, 22™) 
is somewhat restrictive. For instance, it is not obvious how to treat correlators such as 

2iTT 

or even what to expect (as compactification may start to play a role). Indeed, up to a phase, one can write 

(exp( — (/i(x) + h[y) + h{z))) = (cxp( — - h{z) + h{y) - h{z)))) 
and if z' is close to, but at macroscopic distance of, z, we have asymptotically 

(exp( — (/i(x) - h{z) + h{y) - h{z')))) ^\{x- z){x - z'){y - z){y - z')\-^\{x - y){z - z')!' 
Taking z' ^ z yields the heuristic 

(exp( — + /i(2/) + /i(z))) - \{x - z){y ^ z)\-t>\x - y]"! 

which is plainly incorrect (as it is not symmetric in x,y,z). 

We turn to the case where s is not small, and for simplicity we first discuss the case of two marked 
points. Remark that all previous estimates are uniform in s for s in a compact interval of (0, ^). Again we 
use a variational argument; now x,y are fixed and s is varying. Starting from |(;^''(J')~''(^))| = | det(if^K~"'^)|, 
where x = e^*'^* and is implicitly a function of s, we get 



ds 



log|(x''(^)-''(-))| =sRTr(A'^^^) 



and K-y{w,b) = ±2i7rx*^K(w, 6) if 7 crosses the edge (wb). Hence we need to evaluate S-y{w,b) for these 
edges. We have (uniformly in (z, w) in a compact subset of {(z, w) : z x,y;w x,y; z ^ w}) 



1 



S^{z,w) = -Rb 



W — X 



z — yj \w — y J 7r(z — w) 



z-y 



yj 7r{z — w) 



\+o{\x-yr) 
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where ^^73^^ = exp(s'log((u — x)/{u — y))) and u ^ log((u — x)/{u — y)) is a determination branched 
along 7. Since 

z — yj \w — y J z — w z — w \w — x w — yj 
we obtain if 6 w 



S^{h,w)-K-'{h,w) = -Rb e^''^^' -TT^B e-^^(^> = - = +o{l/\x-y\) 

ZTT \ \w — X w — y J J 27r \ \w — x w — y J J 

except if [bw) crosses 7, in which case the LHS is replaced with x^^{S~f{b,w) — K~^{b,'w)). 
Let 7' be a subpath of 7 running from x' to y' , and -£(7') the set of edges of M crossed by 7'. If {vv') = {wh)^ 
(as oriented edges), we have K{w,b) — i{v — v') and K{w,b) = e~"'^"''~"'(''^i(w — v') Thus if Vi,Vi+\ are 
consecutive points on 7' corresponding to the unoriented edge {wb), 

±K{w, b){RB{e'''^'^^a) + RBie-'"^"^^ p)) = ia{vi - Vi+i) - ip{vi - Vi+i) = -ia f dz + ip [ dz 

where ± = + if (bw) G Er{j) and - if (bw) € £^^(7)- Then 

SR y: w)s,{b, n^)=-s^r -L\d^-s^r =1=) d-z + o(i) 

(one may also keep track of the imaginary part, since K{w,b)K^'^ (b,w) ~ p{w,b) is expressed in terms of 
the local gomctry). Wc still need to address the logarithmic singularities at the cndpoints of 7. A discrete 
Cauchy formula together with estimates on Sp{z,w) for z close to x and \w — x\ comparable to \x — y\ leads 
to the estimate 

S,{z, w) - K-U^, = 0(^^^^^ • =0{\x- - 

for |w) — x| < \\x — y\, — x\ < \z — x\ < 2\w — x\. The same estimate holds when interverting x and y. 
Consequently, if we choose 7 such that the lengths of its segments are comparable to their Euclidean length, 
and 7j is the initial segment of 7 (between x and x'), we have 

^ K^{b,w)S^{b,w) = K^{b,w)K-^{b,w) + 0{{\x'-x\/\x-y\)'-'n 

(6«))eE(7i) (bw)eE{ji) 

From the estimate 



-x,xV > ) 2 "\\w-xj n{z-w)J 2 \\w - x J n{z - w) ^' ' ' 



for i Iw — < l-z — x| < — x\, we obtain 



K-U6,»)-K-'(<.,») - ( — ) - S«B (_) +0(|»-.|--) 

for b ^ w (this corresponds to letting y — )■ 00 in the previous estimate). Since the error term is summable, 
it follows that 

Y K^{b,w)K-^^{b,w) = -2s\og\x' -x\+c{x,K,-i,s)+0{\x-x'\^'-^) 
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Moreover, the left handside is unchanged if 7 is another path started at x going through x', as the difference 
can be written as the trace of a commutator. Thus c(a;, A, 7, s) — c{x, A, s). Consequently 

TT{K^S-y) = -2s5R log(y - x)^ + c{x, A, s) + c{y, A, s) + o(l) 

Lemma 23. The kernel (with singularity at x) is local in the sense that K~^(5, w) depends only on A 
restricted to a hall centered at x containing b,w. 

Proof. The discrete exponentials are by construction local. If w is adjacent to x, then K^^{.,w) has an 
integral representation in terms of discrete exponentials (in the variable 6) and consequently K~^{.,w) is 
also local. For a general w, consider a simple path ( ) from X to x' on M\ where x' 

is a face adjacent to w. Let {wibi) = {xiXi+iY. By induction on the length of this path, it is easy to see 
that K~^{.,'w) can be expressed as a linear combination of the K~]j..{-,Wi), which are local. Thus K~^{.,w) 



is itself local. 



□ 



Consider the space {(A,x)} of lozenge tilings rooted at a face (fixing the scale S = 1). As is customary 
for rooted graphs (^3^), we may define a distance by: 

dist{{A, x), (A', x')) = inf{e > : (A, x)iB{x..e-^) isomorphic to (A', x')ib{x' ..e-^)} 

where graph isomorphism are required to preserve roots. We need to take into account embedding data, viz. 
lozenge angles. Thus we may refine the distance as follows: 

dist{{A,x), (A', a;')) = inf{e > : (A, x)|s(a:,e-i) isomorphic to (A', ^-i) and corresponding angles differ by at most e 

Then the set {{A,x) : A satisfies (<!>)} is compact. Indeed, under the condition (4|k), degrees of vertices 
are bounded and the number of vertices in a subgraph is comparable to its area; thus there are finitely 
many isomorphism classes as rooted planar graphs for restrictions {A,x)^B{x,e-'^) (for fixed e > 0). Each 
isomorphism class is parameterised by finitely many angles taking values in compact intervals. 

Using locality of the kernels and their continuous dependence on lozenge angles for a given graph 
type, we conclude that (A, x) i~> c(x,A,s) is a continuous function (for fixed s). In particular it is bounded 
under the condition for A. 

Proposition 24. If s Ci (0, there are constants c(x. A, s) such that 

I (x''(^)-''(^)) I ^ exp(c(x, A, s) + c{y, A, s))\x - y\''^' 
Proof. This is already established for a small sq > 0. We write 

= / diTr{K.yK-^)du 

J Sn 



log 



(^^h{y)-h(x)^ 


-f 




J So 



and we have estimated the integrand. 

Let us now show that c(x, A, s) = c(A, s), which we already know for small s. We need only show 



□ 



(^■<i^h{y)-h{x) I 

if x' £ AI^ is a neighbour of x; in turn, this follows from 

Sp{b,w) - K'^^^{b,w) 



o(l) 
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as y — oo, where (bw) is the edge of M separating x from x' . From Lemma 15 estimates for Sp on the 
macroscopic (|x — scale, and Lemma 16 we get 

Sp{b, w) - K-',{b, w) = 0{\y - 1+^+") 

for rj > 0. This is enough to conclude that c{x, A, s) — c(A, s). 
We now discuss the "general" case: 

n 

{exp{2i7T^Sj{h{yj) - h{xj))) 
i=i 

where si, . . . , s„ G (0, |). The estimates are as in the two point case, so we will simply record the needed 
changes in the computations. This problem is associated to a character p of 7ri(C \ {xi, . . . ,yn}), and the 
kernel Sp inverting K on sections of the associated line bundle. We start from the estimate for Sp (obtained 
from surgery and the two-point case): 

Sp{z,w) = ^Rb {e'''^'"^S{z,w)^ + ^i?^(e-^^('")5(z,u;)) +o(i?-i) 
where all pairwise distances in {xi, . . . , i/„, z, w} are of order R, with 

Siz, w) = n f i^^--^)(^-yA ^ + (^) + oiz - n.) 

n(z — w) -^-^ ^ Vj)\^ ^ ^j) J ^(^ — w) 



SM-^^f[(^,^^^^4^] ' = ^^+fp{w)+0{z-w) 
TT{z~w)l\\iz-yj)iw-x^)J n{z-w) ' 

where the estimate is for z, w away from the singularities, where the Robin kernels are given by 

rpiw) = ^~^j2(— —) 



~ 5 T S j 



Let 7j be a simple path from xj to yj (at macroscopic distance of the other singularities) , and 7^ the subpath 
from x' to w'. Since rJw) ^ - — — near and rJw) ^ — — near -y,, we have 



rp{w)dw = log ((y; - y,){x'^ - Xj)) + 0(1) 



Set 

f rp{w)dw^ lim * rp{w)dw A- —log {{y\ - yj){x\ ~ xA) 

There is no determination issue for the real part of this regularised integral, which is what we need. The 
variational argument above shows that 



9s 1 



^ preg pi^^g 

log I (exp(2z7r Sj{h{yj) — h{xj)))\ — — Sj5R (p nrp{w)dw + SjSR (p Trfp{w)dw 
j=i Jij •'13 

+ dsc{x,A, Sj) + dsc{y,A,Sj) + o(l) 
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and a direct computation yields 



as needed. We have obtained the foUowing result on electric correlators: 

Theorem 25. For weights si,...,S2„ € (— |,^), si + S2 = •■■ = S2n-i + S2n = 0, and singularities 
{zi, . . . , 2;2„} wii/i pairwise distances of order R, we have as R goes to infinity: 



2n \ / 2n 



exp(2z7r^Sj/i(zj)) \ - exp ^c(A, |sj |) J]^ 

\ j = l I \ j = l / l<i<j<2n 



The constants c(A, s) are bounded for A in the class and s in a compact subset of (0, \). Let us 
now discuss the dependence on A (or lack thereof). The locality and continuity arguments for (A, a;) i— > 
c(a;,A, s) = c(A, s) show that, for A satisfying c(A, s) can be approximated within e by looking at any 
ball of A of large enough radius i?(e). 

Consider Kg a regular tiling of the plane by identical lozenges with angles 9, tt — 6. One may glue Ag and 
Agi along a line to obtain a lozenge tiling Ag^gi (as in f2P). Then for small s, 

c{Ae, s) = c{Agfi^,s) = c(Ae/, s) 

by continuity in (A, x). 

More generally, consider the following condition for two rhombus tilings Ai, A2 in the class (<|l): for 
some 6'q > 0, there exists arbitrarily large balls Bi <Z Ai, B2 <Z A2 and a rhombus tiling A3 with angles 
> O'q containing copies of Bi and -B2. This generates an equivalence relation ^ on rhombus tilings; again by 
continuity we have (at least for small s): c(Ai, s) = c(A2, s) if Ai ^ A2. 

The previous argument shows that Ag ^ Agi. As another example, set F to be the triangular lattice 
and At the corresponding rhombus tiling; it contains an hexagonal lattice as a subgraph. Let A' be any 
rhombus tiling obtained from At by star-triangle transformations inside these hexagons. Then At A'; we 
thus obtain an equivalence class with infinitely many elements which are distinct as graphs. 

It is also easy to see that At ^ Ajr/3, since one may glue a half-space of A^ and a half-space space A^/3 
to obtain a planar tiling A. 

The triangular lattice has a two-parameter (up to isometry) family of embeddings with isometric faces, 
leading to a two-parameter family of rhombus tilings. In two steps (as in [H]), one sees that they all belong 
to the same equivalence class. Similarly, the isoradial embeddings of the square lattice (under (♦)), which 
are parameterised by two bi-infinite sequences of angles, fall in the same equivalence class. 

It would be of some independent interest to decide whether ~ is trivial (coarse - and thus c(A, s) = c(s)) 
and if not describe its equivalence classes (building on the structure result of [52]). 



8 Monomers and the Fisher-Stephenson conjecture 

We now consider a planar graph M, derived from a lozenge tiling as illustrated in Figure [T] with "defects" 
consisting in a pair of missing vertices: a black vertex 60 £^nd a white vertex wq. 

Specifically, consider a sequence (S„)„>i of subgraphs of M bounded by a simple cycle on . We 
assume that for all n, 60 and wq are in S„ and that S„ has a perfect matching; and that the inradius of 
say seen from wq, goes to infinity as n — > 00. Moreover we assume that K-i(6,w) ^ K~\6,w) for all 
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b € Mb,w £ Mw (here K„ denotes the restriction of K to S„); this is the case for "Temperleyan" boundary 
conditions (1221 )• It follows that the sequence of the perfect matching measures on S„ converges weakly to 
the measure on matchings on M described earlier. 

We also assume that the limiting measure /i satisfies the following finite energy condition: ^.{{biwi) G 
m, . . . , (6„, Wn) G na} > for any finite subset Wi), 1 < i < n} of edges of M than can be completed in 
a perfect matching of M . 

Following Fisher-Stephenson jl7| . we consider the monomer correlation 

Z(S„ \ {6o, Wo}) 



MonH„(6o, Wo) = 



Z(S„) 



where -Z(S) is the partition function of perfect matchings of the graph S (with edge weight |K(6, w)\ for the 
edge {bw)). 

In general, a Kasteleyn orientation of S„ does not restrict to a Kasteleyn orientation of = S'j = S„ \ 
{&o,wo}. Indeed, if 6o and wq are not neighbours in M, each correspond to a finite face of where the 
clockwise odd condition is violated. In order to obtain a Kasteleyn orientation of S^, one may introduce 
a "defect line" 7, i.e. a simple path 7 on running from a face of M adjacent to 60 to a face adjacent 
to 2/0- Reversing orientations of all edges crossing 7 yields a Kasteleyn orientation of S^. Thus, let us fix 
such a defect line 7 and define K' : ]R*^B\{bo} _^ ^Mw\{wo} K'(w;,6) = -K(w,6) if 7 crosses {bw) and 
K(u',6) otherwise. We denote by K„, KJj the restrictions of K, K' to S„,S^ respectively (throughout, the 
prime recalls the sign change across the defect line). Then [17j : 

^ Z{-n) , det(K;) 
MonH„ 5o, wo = = ± , J " 

Z(s„) dct(K„) 

In order to compare K^, K„, it is rather convenient to extend KJ^ to an operator Mr^ — > M""' by setting 
K'(wo, 60) = 1, with all other new matrices elements K'(i(7o, h), K'(u', 60) vanishing (one may think of an edge 
"handle" connecting directly &o,wo). This does not change the determinant of (up to sign). Then 

MonH„(6o,wo) = ±det(K;K-i) 

Clearly K'^ — K„ has bounded support (vertices adjacent to 7), and does not depend on n (for n large enough) 
on its support. It follows from the assumption of pointwise convergence of K"-'^ that: 

MonM(&o,wo)= lim det(Kj,K-i) = | det(K'K-^)| 

This defines the infinite volume monomer correlation MonM(feoi wo). Under the finite energy assumption, 
MonA/(fcoi Wo) is positive. Let us sketch a direct argument. Consider a simple path from 6 to w on M: 
(60 = ^, Wo, 61, . . . , 6„, w„ = w). Then matchings of containing (wo, ^i),- • ■ j(wri-i, ^n) are in bijection 
with matchings of S„ containing (60, wo), • • ■ , (&n,w„+i) (see [5S] for related arguments), which gives the 
lower bound: 

^/l II. |K(w;o,5i) . . . K(w„_i,&„)| 

MoHAf 6,-u;) > -T-^ — --— pi{ 6o,wo em,..., 6n,w„) e m} 

|K(wo,&o) ■ • ■ K(w„,5„)| 

One may choose the path so that (60, . . . , bn+i) is a simple path on F and is connected to infinity in 
r\{6o,...,6„}. Then 

Ai{(feo,wo) e m, . . . , (&„, w„) e m} = |K(u;o,&o) ■ • • K(w„,6„)(det K"^(fei,Wj))o<ij<„| 

This is positive; indeed, otherwise one could find a non trivial / G (C*^^) vanishing at infinity such that K/ is 
supported on {wo, . . . , w„}, and vanishes at {60, . . . , It is then easy to see that /|r is discrete harmonic 
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except at {bo, . . . , Since the harmonic measure of bn+i seen from infinity in F \ {bo, . . . , 6ri+i} is 

positive, this forces f{bn+i) = 0, then /|r = 0, then / = 0, yielding a contradiction. 

The monomer correlation question bears on the asymptotic behaviour of Mon7v/(6oi ''^o) ^ l^o ~ w^ol ~^ 
eg with Wo fixed. The Fisher-Stephenson conjecture [T7] proposes: 

Monz2(5,w) ~ c\b — w\~i 

In [55], Kenyon analyses the case of a white defect in the bulk and a black defect on the boundary of a 
simply connected portion of the square lattice. This leads to the length exponent for the loop-erased random 
walk. 

More generally, if {fei, . . . , 5„} and {wi, . . . , w„} are n-tuples of black and white vertices of M, one may 
consider the 2rt-point monomer correlation: 

, N Z{En\{bi,...,bn,Wi,...,Wn}) 

MonM(oi,...,6„;wi,...,w„) = lim , 

Let us give a heuristic interpretation of these monomer correlators. If we apply the construction of the height 
function to dimer-monomer configurations on M, with monomers located at the fixed . . . , 6„, wi, . . . , u;„}, 
we obtain an additively multivalued function on M^, which increases by 1 (resp. —1) when cycling counter- 
clockwise around a white (resp. black) vertex. In the absence of defects, the height field converges to a free 
field (in the plane). Interpreting the monomers as discrete versions of magnetic operators for the free field, 
one may expect 

MonM(&i, . . . , 6„; iwi, . . . , w„) - c(: C'i(wi) . . . Oi{wn)0-i{bi) . . . e)_i(6„) :)c 

in some asymptotic regime, for instance as the pairwise distance between insertions goes to infinity (alter- 
natively as the lattice mesh goes to zero). For the planar free field, the magnetic correlator is explicitly 



Zj Z}^ 



j j<k 

and here g = Tr~^, in agreement with eg Corollary |4j This yields the natural extension of the Fisher- 
Stephenson conjecture: 

MonM(^i, . . . , Z2„) ~ cj]^ |zj - Zj\''^^/^ 

where £i = 1 (resp. -1) if Zi is black (resp. white) and J^i — 0- This is supported in particular by Ciucu's 
work (see [5] and references therein). 

The (somewhat stronger) variational form of the statement we shall derive here consists in estimating 
log(MonM(^, w) /MonM(&', w)) within 0{\b— w]^^^^) for 6 black vertices on the same face of M (or similar 
quantities for higher order correlations). 

The relevant family of holomorphic line bundles are the line bundles over E = C \ {zi, . . . , Z2n} with 
holomorphic sections in U of type: 

2n 

su{z) = \{{z-Zi)-^'^g(z) 

4=1 

where g is holomorphic in J7 \ C and vanishing at infinity. 
8.1 Monomer pairs 

In this subsection, we focus on monomer pair correlations, in the framework laid out at the start of the 
section. 
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Let us observe that one may define a dimer measure /ia/' on matchings of M' = M \ {bojWo} as the 
weak hmit of the (weighted) dimer measures on 5^ as n -> oo. It suffices ([23]) to observe that (K^)^^ — 
K~^{K'j^K:^^)^^ , and as n -> oo, converges pointwise to K^^ and (K^K~-^) converges pointwise to an 
invertible (since Mon(6o,wo) > 0), bounded support perturbation of the identity. Let us denote (K')^^ the 
hmiting kernel (or (KJ^^ «;q)~^ need to emphasise the position of monomers). It satisfies K'(K')~^(., w) = 
<5u, for all w G Mw \ {wq}, and (K')-^(6,w) = Oi\b-w\-^). 

Denote by (5o,wi,6i) three consecutive vertices (in counterclockwise order, say) on a face of AI. We 
have: 

Z(S„ \ {bo,wi,bi,wo}) = Z{En \ {6o, wo})|K(u;i, 6i)|" VH„\{6o,t«o}{(^i' "^i) ^ ^} 
= 2(S„ \ {&i,-u;o})|K(u;i,6o)rVH„\{bi,«;o}{(^o, wi) G m} 

where we consider two defects: the monomer wq, and the "trimer" {&o,Wi,5i}. Consequently, 

MonH„ (&i, Wo) 
MonE„{bo,wo) 

and taking the limit as n — >■ cx), 

MonM{bi,wo) ^_^ (K;,^,^^J~^(6i,wi) 
MonM(6o,wo) iK,,wo)~^ibo,wi) 

Thus we simply need to estimate Sbo,WQ — i^'bo wo)~^ precisely enough near the singularity &o- The argument 
is at times somewhat technical; let us briefly sketch the line of reasoning. 

1. Random walk arguments give an estimate on Sp, an inverting kernel for KJ^^ which is bounded (rather 
than vanishing) at infinity. 

2. Bounded functions in the kernel of K'l are classified. 

3. The correct inverting kernel Sbo,wo is constructed, and estimated in the macroscopic scale. 

4. A priori estimates for Sbo,wo when one or both arguments are near singularities are given. 

5. g, an approximation of Sbo^woi-jWi) is constructed, using as building blocks discrete holomorphic 
functions mesoscopically near singularities, and continuous holomorphic functions elsewhere. 

6. 5(61) — Sbg^woi-jWi) is estimated within 0(|6o ~ ''^ol "'^^^) (this requires controlling the leading term, 
which depends only on the local geometry around x, and the first subleading correction which carries 
the "global" information). 



K(wi,6o) 
K(«;i,&i) 



Mh„\{6i,i„o>{(^o, wi) e m} (K' 



Recall that K[,^ is obtained from K by deleting the row and column corresponding to the pair of 
monomers and changing the sign of matrix entries corresponding to edges crossing the defect line 7 which 
runs from x G adjacent to bo to y E adjacent to wq: 

K' . lDA/B\{fto} . mA/»-\{tuo} 

Displacing 7 (with endpoints fixed) results in composing Kbg^wo with diagonal matrices (with ±1 diagonal 
coefficients). Equivalently, one may consider multivalued functions on M corresponding to the character 
/5 : 7ri(C \ {x,y}) {±1} with monodromy —1 around x,y. One may identify elements of M^^f\{bo} -^[fj^i 
functions on AIb vanishing at 60, and elements of M*^»'\{'"o} functions on AIw modulo 6wg. 
We start with by constructing an inverting kernel Sbo,wo for ^ba wq^ least for |6o — wo| large enough. 
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Assume without loss of generality that bo eT^ . We observe that the kernel Sp obtained from the Green 
kernel for the random walk with monodromy at 5o, wq on T vanishes at 60 (see Lemma 17 1 and afterwards). 
Hence ,^,^S'p(., w) = Sw Moreover we have 




Spiz,w) = ^RB I ^/ /^ •"'ii^ f^^\ ]+icc) + oil/\x-y\) 

^ ' " [w — x)[w — y) tt{z — w) I 

for (z, w) in a compact set of {z ^ x,y,w]w ^ x, y, z}. Note that in this case, a closed form expression for 
the continuous Green kernel Gp is easily obtained by going to a double cover of C \ {Bq, wq}. 

The kernel Sp does not vanish at infinity; on the other hand we have some freedom since Sbo,woi-Tw) 
is not required to be holomorphic at wq. In order to correct S to obtain an inverting kernel vanishing at 
infinity, we need to study bounded functions in the kernel of KJ,^ (this is a point where the argument 
differs from the construction for electric correlators). 

We claim that dim{/ : K[,^_^^/ = 0, ||/||oo < 00} < 2. Indeed, K'^^^^^/ = implies ZV/|r is supported 
on the two vertices of T abutting wq, and there are no nonzero bounded p- multivalued harmonic functions 
on r (/|rt is specified by harmonic conjugation and /{bo) = 0). We simply need to construct two linearly 
independent bounded functions in the kernel of K[,^^ . 

We have constructed (Lemma 11 1 g such that K^i yg — outside of wq and 

giz) = RB{{z-y)-'^') + 0{\z-y\-'/') 

We then consider / = glB—Sp{glB) where B — B{y,r), r < ^\x — y\. Plainly, / is bounded and Kh^^^gf — 0. 
Since 



1 / 1 (z — x){z ~ y) dw I z — X I-b(^) 

2j7r Jc(y,r) y/w-y \J {w - x){w - y) z-w ]/ (z - y){y - x) y/z-y 

and the contribution from the conjugate term vanishes, we get 



^(-)=^HV (^~y)(.-^) j-^"^'^-^'"^'^^ 

for z in a compact subset of C\{x,y}). We know a priori that /|p (resp. /|rt) has a limit as z — > 00. Plainly, 
lim^^oo f\r{z) = ~ ^ + o{\x - y\^^^'^) and lim^^oo /|rt(z) = ij^y-x^ o(|x - j/|~^/^). 

We observe that / is also in the kernel of KJ^^ (a real operator); examining the limits of /|r, /|rt at 
infinity shows that /, / are linearly independent. Moreover, there is a unique kernel Sha.wa inverting K[,^ 
and vanishing at infinity. By uniqueness, this kernel is real. It may be realised as 



560,^0(6, w) = Sp{b,w) - 5R(a(w)/(6)) 

(w—x) (w — y) 



where a{'w) is fixed by the behaviour as 6 — > 00: a{w) = — ^/ (y^_^)(^_y) + o{l/\x — y\^^'^)- Consequently, 
we obtain the estimate 




/ N -r. 1 iz-x)(w-y) e*'^^"') \ , , , „ 
S,o,^oiz, -) = t^Rb I \l (,_^)(^_,) • J + (-) + o{l/\x ~ y\) 

for (z, w) in a compact set of {z ^ x,y,w;w x, y, z}. 

Let us denote by {bQ,wi,bi,W2) the vertices of the face x of M, listed counterclockwise. We want to 
estimate Sbo,wo{wi,bi) within 0{\x — j/|^^^^). We shall need basic discrete holomorphic functions with 
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monodromy x — around the singularity x G M'''. Set g = K^"^(., wi), /-i = at x = ^1; h = fi-i/r 
where t = (61 — x)/\bi — x\. We Hst some useful evaluations (Lemmas [lO 11 ). 

r(i) / e'''('"i) 



2 \ n[b — w) y wi — X I 



2 sin ^0 



5(60) 

2 sm t^i 

hib) = V2Ti~^)RBiVb^) + Oi\b ~ 

h{bo) = -t7r{2ibo - x)y/^ 
- -7r(2(6i - :r))i/2 

where 9i — arg(^5^)i and we assume that the defect line 7 does not cross (boWi), (biwi). 

Let us review some a priori estimates on Sbg^waib^w) for b close to cc, and w close to x^y,oo respectively. 

Lemma 26. There is e > such that: 

1. IfK^i,J = OmB{x,r), /(60) = 0, then f{b) = 0{i\b - x\/ry\\faBi^,r)\\oo) forbeB{x,r). 

2. There is a unique kernel Sbg s.t. Sb„{-,w) has monodromy (—1) around x, K_iS'b„(.,w) = 5w, and 
Sbo{-,w) vanishes at bo and infinity. 

3- Sbg.wQ {b, w) — 0{\w — xj^"'"^^) if b adjacent to x, \w — x\ < 5 |a; — y\ 

4- Sbo^wo{b,w) = 0{\w - y\^/'^^^\x - yl"^/^) ifb adjacent to x, \w ~ y\ < l\x - y\ 

5- Sbg.wo {b,w) = 0{\w — x\^^\x — y\^'^) if b adjacent to x, \w ~ x\ > ^\x — y\,\w — y\ > ^\x — y\. 

Proof. 1. This follows from observing that /|r is harmonic, with monodromy (—1) around x, and writing 
/|pt as the harmonic conjugate of /|p vanishing at fog. 

2. Uniqueness follows from the fact that the space of (— l)-multivalued bounded holomorphic functions 
is spanned by a single function f_i which does not vanish at bg. By uniqueness this kernel is real. As 
before (when wq is at finite distance), one can write Sb„{.,w)\i' as a first difference of the Green kernel 
on r (with monodromy —1 around x) and extend it to by harmonic conjugation. This leads to the 
a priori estimate 

Sb„{b,w)^Oi\w\-\{\b\/\w\rA{\w\/\b\Y)) 
eg if \b — w\ > ^\w\ (setting a; = for notational simplicity). 

3. Let us write 

Sbo,wo{-,w) = lBSbo{.,w) - Sbo.w„{K„,woi'^BSbo{-,w))) 

where B = B{x, l\x - y\). We have K'^^^^ (IsS't,, (., w))(w') = 0{\w\~^{\w\/ Rf) for w' G dB (where 
R=\x-y\), &nd Sb,^^,{b,w') ^ 0{R-^'-') by 1. Thus Sb,,^,{b,w) ^ 0{\w\-^-'). 

4. Here the reference kernel Sy is constructed from the Green kernel which has monodromy (—1) at y and 
is such that Sy{b,w) = 0{\w — y\~^) ii\b — y\ and |6— w| are comparable to jw — t/|. We can construct 
two functions fy, fy with monodromy (—1) at y which are holomorphic except possibly at wq and 
with expansion fy = _Rb((6 — wq)^^^^) + 0{{b — wq)^"^^^). Correcting Sy{.,w) by an appropriate linear 
combination of fy, fy, we obtain a kernel such that Sy{b, w) = 0{\w — y\^^{\b — y\/\w — y\)^"^^^'^^) 
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for \b — y\ > 2\w — y\, say. 

Then we may represent Sb„.wo( ■■'"'') 

Sbo,wo{;W) = lBSy{.,w) - Sbo,wo{'^BSy{.,w)) 

with B = B{y, — y\) to obtain for b close to x: 

Sbo,v.o{b, ^v) = 0{\x - yr\w - y\-\\w - y\/\x - y\f/'-') 

5. We simply write 

Sbo,woi-,w) = IsK {.,W) - Sbo,wo{^B\<r^{;W)) 

with B = B{x, l\x - y\) U B{y, \\x - y\). 

□ 

We now want to construct an approximation of Sbo,wo{-> ^i)- We are looking for such a function g of the 
form: ^ is a linear combination of g, g, h, h in B{x, r); 



g{b) = ^Rb 



' [b - kt')[tri - wq) 
7r(6 - wi) y {wi - bo){b - wq) 



in {B{x, r)L)B{y, f)y. In B{y, f) we set g = ^{"ffy) where fy is a function such that fy{b) = RB{{b—y) ^/^) + 
0{\b — and fy is (— l)-multivalued around y and is holomorphic except at wo- Set R = \x — y\. 

Wc need: .g(5o) = 0, (Kj,^^ ^jj,g)(wi) = 1, and the values of 5 across dB(x.r) differ by o{^/R): this 
defines uniquely g as a linear combination of h, h and fixes fi. Since g satisfies the same conditions, we 
may set g = 5R(ag + in B{x,r). Wc get fi = fe"'(-'"^'>r{l/2),J^^^. 

Then g — Sbo,wo O^bo wo9 ~ ^wi) = Sbf^^woi-) '^i)- We need to fit a, j3 in order to make the error term small 
near 60. Since 



{b - bo){w-i_ - Wq) b-bo (11 1 , ^/li, l/IJ, |2- 



h V( ^ \/ ^ \^ +0 ,+0{\b-w^\l\b-WQ\') 

b — wiy [wi — bo){b — Wq) \ Wi — bo \b — Wi 2 wq — b 

^ {l + 0{\b-wo\/R)) 



^y{wi - bo)ib-wo) 
we require: 5R(a) = 1, 5R(ag((6o) + l3h{bo)) = and 



2 2^^^/Wl — x 2{wo — bo) 



or /3 = ae"^("'iV(87r v/2(w;i - x){bo - wq)) = az with z = 0(l/i?). Let us remark that = | does 

not depend on the local geometry (while ^{g{bi)) docs). From di{ag{bo) + (3h{bo)) = 0, we get 



^{g{bo) + mbo)) = \^{a) 



and thus a = 1 + 2i^{g{bo)) = 2ig{bo) + 0{1/R), and by iteration: 



a = 2ig{bo) + 4m{ig{bo)zh{bo)) + 0(1/ R^ 

and we have 



g{bi) = 3fJ(a.9(6i) + ^h{b,)) = 2^{ig(bo)g(bi)) + m{ig{bo)zh{bo)){--) + '•R{2ig{bo)zh{bi)) + 0{\/R^) 
= 2^{igibo)g{bi)) + ^{2ig{bo)z{h{b^) - h{bo)) 
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Then 0{r-'^/^ / R + r'^^'^ +r^/^ /R'^) = 0{^/¥/R + r~^/^) on dB{x,r). This contributes to a bound 

of order 0{r^^^^^'^ R^^ + r^'^l"^^^ + r'^l'^^^Rr'^) on the estimation error (5(6) — S\,^^w^{h,W\)) for 6 within 
0(l)ofx. 

In R)\B{x, r), we have K'j,^ ^^^ g = 0{\b~x\^^^'^^^), contributing to an error of order 0(X]feLr k^^^^^^^'^) — 
(9(7--i/2-2-£-)^ In B{y,R) \ B{y,f), we have (KJ,^ u;o5)('"') ~ 0{\w ~ wo\^^^'^^^), contributing to an error 
of order 0{J2k=f k-^^^'^k^/^-'' R-^/^) = 0(f-i-^i?-3/2) Qn dB{y,f), we estimate K^^^J = 0(f-3/2 _^ 
leading to an error ©((f"^/^ + f'^/'^ / R^^/"^^^ R^^/"^). Finahy, outside of B{x, R) U B{y, R) we have 
iKo.w.M'^) ^ 0{\w\-'^^/R) and an associated error of order 0{YX^Rk^^VR-R^^) = OiR^/"^-^-^). 
Summing up, by setting r = R ' , r = o we obtain: 

5'fco,»o(6i,i(;i) = 5(61) + O{\bo - woT^""). 



We are interested in exchanging the roles of bQ,bi. We have 

^60.-0(^1: ^1) = 23?(i^g(6i)) + 3? (2i^zihibi) - /i(6o))) + 0{R-'~') 



5, 



Observe that 23^(45(60)5(61))) = .9(^1) — 5(^0)1 which is real (since 3(.g(6o)) = 5(5(61)) = |), and that 
h(bi) ~ h(bo) = -27r(2(6i - x))^/^. An examination of the local geometry yields ie*''("'i) ^'^M^M^ ^ 



-^(61 — 60). Consequently, 

1 + •R{2iz{h{bi) - h{bo)) + O{\bo - wq 



= i-l^(^^]+oi\bo-wor-n 



Ion A/ (61, Wq) _ Sbo,woibi,Wi) _ ^ ^ sa(^A„ihih \ hih \\ , „.. 

^Or\M{bo,Wo) Sb^,wa[bl3,Wl) 

= 1 - , 

2 V ^0 - Wo 

and thus we have obtained: 

Proposition 27. There is c(A, w) > such that 

MonM{b,w) c(A,w;)|6- wp^/^ 

as b ^ 00. 



Moreover c is a continuous function of the face-rooted rhombus tiling (A,w), and is bounded away from 
and CO on (4|k). Presumably, c depends only on A (and thus is constant on a large class of rhombus tiling, 
see the discussion after Theorem 25). To see this, it is enough to check that 

Mon(6, w') ^ 
Mon(6, w) 

as 6 — >■ 00, if w, 6', w' are consecutive vertices on the boundary of a face of M. In turn this is equivalent to 

Sb,wib',w') 
Sb,w'{b',w) 

as 6 00. This seems to require another ad hoc local computation, which we leave to the dedicated reader. 

Notice however that if M is the square (or rectangular) lattice, one can interchange the roles of black 
and white vertices and consequently c(Z^, w) — c{I?). 
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8.2 General monomer correlators 



In this subsection we indicate how to extend the previous argument to the general case of 2m monomers 
61, . . . , 6m, wi, . . . , Wm in the plane, which does not involve any substantial additional difficulty. We denote 
b = {61,..., 6™}, {wi,...,Wm}- 

Let Xi (resp. yi) be the face of M adjacent to bi £ Mb (resp. Wi G Mw)- Let 7^ be a defect line 
running from Xi to yi on M^; we assume that the 7i's are disjoint and that the pairwise distances between 
singularities are of order i? 3> 1 . Let 

J^' . ][jA/B\{f)i,...,6„} ^ jjA/vrXfii'i,---,-!"™} 



be the operator obtained from K' by removing rows and columns corresponding to monomers and changing 
signs of entries corresponding to edges of M crossing one of the defect lines. Equivalently, one may think of 
K' as mapping {/ e : /(61) = • • • = /(5„) = 0} to M*^^"' /]r{»i,-.«'™}. 

By surgery (Lemma [s]) and induction on m (as we did for multiple electric correlators), we get that for 
R large enough, there exists a unique inverting kernel Sb,w for K' vanishing at infinity and that 

Sb,^{z,w) = ^B.B [e'''^^^S{z,w)') + (cc) +o(l) 
where all pairwise distances are of order R and S is the continuous kernel 



(z - Xt){w ~ y^) 



s{z,w) = n , . w . , , . 

^J^\[z - yi)[w - x^) J n(z-w) 



1/2 



The corresponding Robin kernel is 

r{w) = lim S{z,w) — 

so that 

1 



1 



tt{z — w) 



9-7r ^ — ' \ HI — 



Stt \w — Xi w — yi 



f{xk) = lim r{w) 



-y.(— 



2TT{'W~Xk)J 2TT{yk~Xk) 2'K^^\Xk-Xi Xk-yi 



If (60, w, b'o) are on the boundary of the face Xk (with 60 = bk), a straightforward extension of the two-point 
argument yields 



Sb,n,{b'o.w) = 2'R{ig{bo)g%)) + ^ 



2^2{w - x) 

= 9{bo) - gibo) + {g{bo){bo - 6;,)^f(a;,)) + OiR-'-^ 



fixk) +OiR-'-') 



(compare also with Lemma 20). Let b' — (61, . . . ,&n) with &q substituted for bk = bo. Then taking into 
account g{bi) — g{bo) S K, we get 

±%4^ = 1 + 5? m - bo)nfixk)) + 0{R-'-n 



From 



Sb,w{b'o,w) 

(: Oi{w,) . . . Oi{w,rr)0-i{b,)0-i{b^) :)c = c 



I\i<Abi - bj)iwi - Wj) 



1/2 
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we get the following variation when bk = 60 is replaced with b^: 

Consequently, we have: 

Theorem 28. Let 61, . . . ,6^? wi, . . . ,w,n be m.arked black (resp. white) vertices on M, with pairwise dis- 
tances of order 1. Then there is Cm{^,'w) > such that 



Mopm^w) = c„(A,«)) 



1/2 



(1 + 0(1)) 



Again in the case of the square lattice, one may switch colours, so that Cm{^,w) = Cto(Z^)- Presumably 
Cm(Z^) = ci(Z^)™, though that requires an additional argument. 

Let us point out a by-product of the method. In the course of the argument, we have considered trimers 

of type {60, Wo, &o}, three consecutive vertices on the boundary of a face / of M. Then we have for example 
the following expression for the trimcr-mononicr correlator: 

MonM{{bo,wo,bQ},wi) = MonM(bo:Wi)\Sb„^wAbo,'Wo)\ = MonM{bo,wi)\Sba,wi{bo,wo)\ 

and |S'(,(,,^i(&o,U'o)| = 15(^*0 - ff(^o)|(l + o(l)) as \bo - Wi\ ->■ 00, where 5(60) - g(bo) depends only on the 
geometry of the face /: 

9%) - 9{bo) = 2 cot arg (^^-^ j - - cot arg (^^-^ j / 
where x = ig center of /. In the general case, we have 

'cot(6ii) +tan(6'i) 



MonM{{bi,w[,b[}, w^, 6^}, Wi, . . .,Wm) ~ Cm{^,w) 



Ui<j{bi - bj){wi - wj) 



Uijtjibi - Wj) 



1/2 



where e (0, |) is an angle of the right-angled triangle {6j,«;^,6^}. 

Let us phrase a more general conjecture, following Ciucu. For simplicity we consider the case Af = Z^. 
An islet / is a finite subset of vertices of bounded by a simple loop on (Z^)'''. Its charge is e(/) = 
|/ n Mw\ — |/ n Mb\- The conjecture is that an islet / is a discrete version of a magnetic operator Oe(/) in 
the sense that 



(1 + 0(1)) 



MonM{h+Xi,...,In + Xm) C{h,...,ln) \Xi - Xjl'^^^^^'^^^^^'^ 

l<i<j<n 

as the pairwise distances \xi — Xj\ go to infinity (we require here X)je(/i) = 0). 

Lastly let us point out that the surgery argument enables to analyse at essentially no additional cost 
mixed magnetic-electric correlators, the simplest of which is 

{Oi{w)0.i{b) exp(2i7rs(</.(/') - </.(/)))) 

where s E (0, ^) and the pairwise distances between insertions b,w,f',f go to infinity (and Oi{w), 0-i{b) 
represent monomer defects). Of some interest are the coincident magnetic-electric operators. At the lattice 
level, one may consider the above correlator with w on the boundary of / and b on the boundary of /'. For 
general s, their analysis seems to require additional arguments. Note however that in the coincident case, 
for s = i, 

(Oi(w;)0_i(6) exp(z7r(</>(/') - ,^(/)))) = ±K-\b,w) 
the asymptotics of which underpin the analysis of all other correlators. 
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